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ABSTRACT 

We explore the use of simple star-formation rate (SFR) indicators (such as may be used in 
high-redshift galaxy surveys) in the local Universe using [O II], Ha, and u-band luminosities 
from the deeper 275 deg 2 Stripe 82 subsample of the Sloan Digital Sky Survey (SDSS) cou- 
pled with UV data from the Galaxy Evolution EXplorer satellite (GALEX). We examine the 
consistency of such methods using the star-formation rate density (SFRD) as a function of 
stellar mass in this local volume, and quantify the accuracy of corrections for dust and metal- 
licity on the various indicators. Rest-frame w-band promises to be a particularly good SFR 
estimator for high redshift studies since it does not require a particularly large or sensitive ex- 
tinction correction, yet yields results broadly consistent with more observationally expensive 
methods. We suggest that the [O II]-derived SFR, commonly used at higher redshifts (z~l), 
can be used to reliably estimate SFRs for ensembles of galaxies, but for high mass galaxies 
(log(Af*/M0)j>lO), a larger correction than is typically used is required to compensate for 
the effects of metallicity dependence and dust extinction. We provide a new empirical mass- 
dependent correction for the [Oil] -SFR. Note: This astro-ph version has been updated 
to correct the typographical errors in equations 2, 7, and 9 of the originally published 
paper, as described in the MNRAS Erratum. 

Key words: galaxies: evolution - galaxies: general 



1 INTRODUCTION 

The most useful observables in constructing a picture of galaxy for- 
mation are likely to be: a galaxy's star formation rate (SFR), its stel- 
lar mass, and the epoch at which it is obse rved. Stellar masses are 
now routinely measured out to z~4 (e.g., Marc hesini et alj[2009h 
by fitting spectral energy distributions (SEDs) of stellar population 
synthesis models to multicolour broadband photometry. Although 
the uncertainties in stellar mass estimates are dominated by evo - 
lutionary uncertainties within the models dMarchesini et alj|2009l) . 
most workers use the same models and similar passbands for the 
photometry (or the differences between models are at least well un- 
derstood) and thus the comparison between stellar masses used in 
different works is relatively straightforward. The situation is some- 
what less uniform for SFR measurements. Indeed, most of the sys- 
tematic errors affecting SFR measurements are likely to be func- 
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tions of stellar mass | Brinchmann et al. 2004|; iKewlev et al]|2004l : 
ICowie & Barge3l2008l : IPannella et al]|2009h . 

Over the last decade, measuring the star formation rate density 
(SFRD) of large samples of galaxies out to z~6 has become rou- 
tine, and dozens of such measurements now exist ( Lilly et al. 1996; 
i Madau et al]|l99fj|; ISteidel et al.|[l999l ; iHopkins & BeacomlkoOr: ; 
iReddv & St eidel 2009). In this time of 'survey science', statistical 
samples are commonplace and the largest remaining uncertainties 
in such studies are systematic in nature, for example, the choice of 
SFR indicator. Common proxies for SFR include the luminosity of 
emission lines such as Ha, or the UV continuum, which are sen- 
sitive to the ionising flux from young stars; or mid- or far-infrared 
radiation (MIR, FIR) which measure the fraction of this flux ab- 
sorbed and re-radiated by dust. 

Each SFR indicator possesses its own strengths and disadvan- 
tages. The luminosity of the Ha recombination line is relatively 
simply and directly coupled to the incident number of Lyman con- 
tinuum photons produced by young stars, and hence is propor- 
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tional to the SFR dKennicuttl 1 19981) . although there is some de - 
pendence on the metallicity of the gas tCharlot & Long hetti 2001). 
Provided the line can be adequately corrected for stellar absorption, 
the largest uncertainty is the correction for the presence of dust. At 
higher redshifts, z > 0.4, the Ha line passes out of the optical 
window and the most commonly-used emission line in the range 
0.4 < z < 1.5 has been [O IDA3727. [Oil] luminosity is more 
strongly affected by dust than Ha. Moreover, this collisionally- 
excited line is very sensitive to metallicity. Since the emission lines 
necessary to accurately measure dust extinction and metallicity are 
not available when [O II] is used (otherwise the Ha line would be 
used in preference to [O II]), alternate approaches to empirically 
correc t the [O II]-SFR for these dependencies have been adopted 
(e.g.. Moustaka s et alj2006l : [Weiner et alj2007h . Another observa- 
tional limitation of indicators based on spectroscopy, particularly 
for low redshift, is the need for an aperture correction. The spec- 
troscopic aperture may only sample the innermost regions of the 
galaxy and thus an extrapolation is necessary to sample the total 
light from the galaxy. This may be further complicated by popula- 
tion gradients across the galaxy and the fixed angular size aperture 
sampling a different fraction of a given galaxy at different redshifts. 
A compreh ensive examination of ap erture effects within SDSS data 
is given in Brinc hmann et al. 

The rest-frame UV luminosity is straightforward to measure 
in high redshift surveys. Indeed, it is normally produced as a by- 
product of the survey data, since all that is required is a relatively 
blue optical band (depending on the redshift of the source). The 
UV luminosity provides a measurement of the continuum radi- 
ation from young stars, again after it has been extinguished by 
dust. Although se nsitive to dust, the UV is relatively insensitive 
to metallicity (e.g.,[G]azebrook e t aljl999h . Additionally, a correc- 
tion must be made for the contribution to the UV luminosity from 
more evolved stellar populations. This correction becomes larger 
towards lower redshifts, as the Universe ages and so do the typi- 
cal ages of the galaxies considered. The extinction due to dust may 
be estimated from the slope of the UV continuum (which is easily 
measured from two or more filter photometry). This method was 
originally applied, with cons iderable uncertainty, only to starburst 
galaxies (Meur eret al.ll 1999h . but recent work has found c ompa- 
rable relations for normal galaxies (e.g. JCortese et al1l2006b . The 
effect of dust on these indicators can be large and varies consider- 
ably from galaxy-to-galaxy. In the local Universe, the dust extinc- 
tion in Ha can be -2 mag and in the UV 0-4 mag dKennicuttl 19981 ; 
ICalzetti et alj200Cb . This scatter may be partly driven by galaxy in- 
clination and dust geometry effects, meaning that even for galaxies 
with the same dust content, a different effective attenuation may be 
observed. 

Furthermore, the different indicators are sensitive to the pres- 
ence of young stars for different lengths of time: the ionising ra- 
diation responsible for the Ha line measures SFRs on timescales 
of ~10 Myr, whereas UV-bright stars which are considered in the 
SFR derived from UV luminosity have lifetimes ~100 Myr. Simi- 
larly, the sensitivity of the ultraviolet continuum and emission line 
luminosities to stars of different masses (J> 5 and J> 10 Mq re- 
spectively) means that the derived SFRs are sensitive to the IMF. 
This also means that if the assumption of a universal IMF is in- 
correc t, this will lead to differences between these SFR indicators 
(e.g jMeurer et alj|2009l) . 

Thus, the SFR is difficult to estimate accurately for an indi- 
vidual galaxy. One solution is to measure the SFRs for large en- 
sembles of galaxies and compare their average properties in or- 
der to smooth over some of these effects. For example, comparing 



globally-averaged SFR quantities estimated from different indica- 
tors with different sensitivities to dust, metallicity, etc. should show 
the same results, on average, if the corrections for these effects 
are broadly consistent. Now, the potential systematic error of most 
concern is the timescale to which the different indicators are sensi- 
tive. However, provided the timescales of all the indicators used are 
short compared with the timescale over which the cosmic average 
SFR is changing, this difference can be negated. 

Since most of the variables affecting the estimate of SFR cor- 
relate with mass, e.g., metall icity via the mass-metail icity relation 
(Tre monti et al] |2004l) . dust dBrinchmann et alj|2004l) ; an obvious 
first step is to determine mass-dependent corrections to the SFR 
indicators. Furthermore, there has been great activity recently in 
measuring how SFR (or some related quantity such as SFRD or 
SFR/stellar mass) varies as as a function of stellar mass, and how 
this evolves with redshift. Most of these results have been consis- 
tent with a picture in which the termination of star-formation pro- 
gresses from higher mass to lower mass galaxies as t he Universe 
ages, so-called cosmic downsizing dCowie et al.lll99ot) . Obviously 
any residual systematic errors as a function of mass on the SFR 
indicators used could pose a serious problem when studying such 
relatively subtle mass-dependent effects. 

To date, most high-redshift studies have used locally- 
calibrated tracers of SFR, under the assumption that the calibra- 
tion is still valid at higher redshift. Most local SFR calibrations use 
a simple conversion between luminosity of some SFR proxy and 
SFR and assume that this is valid for all masses. However, as we 
will show in this paper, a better approach is to allow for a mass- 
dependent conversion factor which accounts for differences such 
as the average dust extinction as a function of galactic mass. Dust 
obscuration at higher r edshift is an open que stion . For example, in 
two re cent z~2 studies. fPannella et al. (2009) and Redd v & Steidel 
d2009l) both find evidence for mass-dependent extinction, but the 
former find an extremely large dust correction is necessary (much 
larger than that measured locally, § 16,2,41 . whereas the latter find 
that the amount of dust extinction required is a decreasing func- 
tion of redshift. P art of this difference ma y be due to sample se- 
lection in that the iReddv & Steidel d2009l) sample selects Lyman 
Break Galaxies, which may favour systems with lower UV extinc- 
tion, but this still goes to show that the evolut ion of dust proper- 
ties is not well constrained . Other works (e.g. jMartin et al.ll2007l : 
Igl esias-Paramo et al l2007h at z~ 1 find similar evidence for mass- 
dependent extinction using the ratio of IR to FUV luminosity to 
estimate extinction. 

Regardless, a sensible null hypothesis would be to assume 
that the higher redshift relations (such as the mass dependence of 
extinction) follow the local relations, and derive SFR calibrations 
and corrections based on local data. Using a local mass-dependent 
correction at higher redshift is clearly preferable to ignoring the 
mass dependence of SFR indicators altogether. As a first step in 
understanding how the various SFR indicators commonly used at 
higher redshift relate to one another, we have undertaken a com- 
parison between these indic ators in the loca l Universe. The Sloan 
Digital Sky Survey (SDSS. lYork et alj|2000t) provides high qual- 
ity optical photometry and spectra for a well-understood sample of 
galaxies. In addition, several groups have already performed de- 
tailed measurements of SFRs using state-of-the-art techniques (fit- 
ting spectrophotometric data to the latest stellar evolutionary syn- 
thesis and photoionisation models) which may be used as a ref- 
erence against which simpler methods, such as those necessarily 
employed at higher redshift, may be compared. 

The main aim of this paper is to examine the consistency be- 
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tween different star-formation indicators by exploring the SFRD 
estimates obtained with the different tracers and to look for possi- 
ble systematic differences as a function of mass. We examine these 
in the local Universe, but end with a brief exploration of how these 
may affect higher redshift surveys. This paper is designed to serve 
as a local reference, both for understanding how the various es- 
timators compare with each other and to examine evolution when 
used in conjunction with higher redshift surveys. We proceed in the 
following way. In §2 we describe the SDSS optical imaging and 
spectroscopy and GALEX UV photometry. §3 describes the SFR 
indicators examined. We use the Ha and [O II] emission lines and 
UV photometry from the SDSS u-band and GALEX FUV. For each 
of these indicators, we explore a range of different assumptions 
for the dust correction, e.g., a nominal constant extinction, extinc- 
tion estimated from optical emission lines (or, for example, the UV 
slope). We present results for our 'best estimator' for each indica- 
tor throughout the paper, but also show how these estimates change 
under the various different assumptions commonly made. Many of 
the prescriptions we use are well known in the literature (particu- 
larly the UV), but we will show that current [O II] indicators are 
far from optimal and first we will need to correct this. In §4 we 
present our new empirical (mass-dependent) correction to [Oil]; 
§5 examines the emission line luminosity functions and §6 presents 
the the SFRD results, convergence tests for the limiting SFR nec- 
essary to accurately estimate the SFRD, and also briefly considers 
the application of these to a higher redshift survey. A more detailed 
z~l analysis drawin g on these local results will be presented in 
iGilbanketafl feoiO). In 87 we summarise our conclusions. 

All magnitudes are quoted on the AB system unless otherwise 
stated, and we assume a cosmology (h, Qm, Q\) = (0.7, .3, 0.7). 
Throu ghout, we convert all quantities to those using a Kroupa 
(2001) universal initial mass function (IMF), except for \ 



2 DATA AND METHOD 

The SDSS is a project dYork et al.ll2000h that has imaged 10 4 deg : 
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and obtained spectra of 10° objects dAbazaiian et al]|2009l) . The 
imaging was obtained using five broadband filters by drift scan- 
ning with a mosaic CCD camera on a 2.5-m telescope dGunn et all 
Il998h . Various samples were selected for spectroscopic targeting 
using a 640-fibre fed spectrograph on the same telesco pe. The sam- 
ples i nclude an r < 17.8 main galaxy sample (MGS; lstrauss et all 
2002) and a colour- se lected r < 19.5 luminous red galaxy (LRG: 
Eisenstein e t al. 2001) sample. 



The imaging survey covered ~ 8000 deg of the NGR Only 
three stripes (scans each 2.5-deg wide) were observed in the SGP 
as part of the 'legacy survey'. One of these, known as Stripe 82, 
here defined as RA from —50.8 to 58.6 degrees and DEC between 
— 1.26 and 1.26 deg rees covering an area of 275 de g 2 , has been im- 
aged multiple times ( Adelman-McCart hv et alj2007l) . In this paper, 
we only use data from the deeper Stripe 82 region in order to probe 
to the lowest possible stellar masses. For our imaging catalogue, 
we use the coadded -catalogue technique and data set described in 
iBaldrv et alj d2005t) . At the time, there were between 6 and 1 8 re- 
peat images depending on the sky position. The data set uses av- 
erage fluxes determined at the catalogue level, providing higher 
S/N magnitudes and colours. This is sufficient for our purposes 
because the surface brightness limit of SDSS single-scan imaging 
dBlantonetafll2005l) is probably not missing any sig nificant de- 
tectio ns for galaxies with stellar masses > 10 s 5 Mm dBaldry et al.1 
2008); and while there is a catalogue derived from Stripe 82 coad- 



ded imaging in DR7, the reliability of the earlier coadd at the cata- 
logue level has been thoroughly checked. 

The catalogue-level coadd is described in § 3.1 of lBaldrv et al.l 

d2005h - This was for the u-band Galaxy Survey. The method ap- 
plied here is the same, except the requirement for fluxes with S/N 
> 3 is applied to the g- and r-band Petrosian measurements but not 
the u-band. For this paper, we select sources with r ps f — r mo£ j e i > 
0.2 and rp e t ro < 19.5. In addition, objects were selected to be 
not saturated, which is a standard flag check to remove artefacts 
and sources near bright stars. This produced a catalogue of 204 093 
objects. 

In addition to the multiple-scan imaging, various additional 
spectroscopic targeting programs were implemented on Stripe 82. 
In addition to the MGS and LRG selections, galaxy targeting in- 
cluded the it-band Galaxy Survey, a low-z (< 0.15) and high- 
2 (> 0.3) selection, and fainter LRG selection. These are de- 
scribe d in § 3 of the SDSS DR4 paper dAdelman-McCarthv et al.l 
2006). Notably the low-z program also includes random sampling 
to r < 19.5. From the SDSS DR7 data release, 73 035 spectra 
were matched to the r < 19.5 catalogue used in this paper. Of 
these, 46 279 are for sources fainter than r > 17.8. While these are 
fainter targets than the MGS, the redshift success rate is still high 
and we do not make any correction for missed redshifts. 

The targeting completeness of our sample is nearly 100% at 
r < 17.8 and drops significantly above that. To model this, given 
the heterogeneous selection, we divide the sample by r-band mag- 
nitude and by the colours u — g and g — r similarly to lBaldrv et al.l 
(2005). Each bin contains 60-166 objects and the completeness, d, 
of sources within a bin is defined as the number of objects with an 
SDSS spectrum divided by the number in the bin. Colour-colour 
completeness maps are shown in Fig. [T] Where the data are only 
available for spectroscopy in the DR4 release, the Ci values are re- 
computed accordingly. 

We select SDSS galaxies in the redshift range 0.032 < z < 
0.2, where the low redshift cut is imposed to ensure that [O II] falls 
within the spectral range of the SDSS spectrograph. The MPA 
group has produced value-added data products for the SDSS, how- 
ever not all products are available for all generations of the survey. 
Currently, emission line fluxes are available for DR7, but stellar 
masses and spectral classes (described below) are only available up 
to DR4. As such, we choose to derive the missing properties us- 
ing simple scalings (which we detail below) from the existing data 
where necessary. The Stripe 82 data are mostly covered by DR4; 
out of the 48 211 galaxies within our redshift range, 43 155 have 
MPA-derived quantities from DR4. 



2.1 Stellar masses 

Stellar masses, M« , for D R4 galaxie|3are measured as described in 
iKauffmann et alj d2003bl) . Briefly, these use the z-band photometry 
to measure the galaxy's luminosity. A galaxy's mass-to-light ratio, 
M/L, is highly dependent on its star-formation history, SFH. Con- 
straints on the SFH are obtained from the galaxy's location in the 
D n (4000) - H8a plane, using a sui te of stellar population synthe- 
sis models dBruzual & Charlotl2003l) to marginalise over a range of 
star-formation histories. 

For the galaxies in DR7 not present in DR4, these stellar 
masses are not available, and so we use a lumin osity and colour- 
based estimate of the mass-to-light ratio (see e.g.. lBell et alj|2003h 



1 http://www.mpa-garching.mpg.de/SDSS/DR4/ 
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Figure 1. Completeness maps: average targeting completeness as a function of u — g and g — r for five magnitude ranges shown by the coloured squares. The 
observed colour-colour bivariate distributions are represented by the black contours. The completeness patterns in the lower panels are primarily the result of 
the u-band, low-z and high-z galaxy selections (annotated in the lower-left panel as a guide). 



fitted empirically against the Kauffm ann et alj d2003bh masses for 
this small subsample of objects. The ste llar mass can be e s timated 
to within 0.21 dex (la scatter) of the iKauffmann etail d2003bl) 
mass using 

log(M*/M ) = 0.480(g - r)o - 0.475M. - 0.08 (1) 

where (g — r)o is the rest-frame (g — r) coloufl and M z is 
the rest frame absolute z-band magnitude. Rest-frame magnitudes 
are estimated by fc-correcting observer-frame magnitudes using k- 
corrections determined u sing the public code KCORRECT v4_l_4 
(Blanton & Roweis 2007). Other colour combinations give slightly 
larger scatter^ 

We also estimate stellar masses for all galaxies using SED fits 
to the ugriz photometry at the spec troscopic redshift using the 
technique of Glazebrook et al ] < l2004l) . such as might be done at 
high redshift. Our results are negligibly affected if we substitute ei- 
ther th e simple colour masses or the SED-fit stellar masses in place 
of the Kauffma nrTetai] d20 03b) mas ses. So, throughout this paper 
we use the IKauffmann et al.l d2003bl) masses, with colour masses 
used for the small fraction not possessing measurements of these 
(DR7 galaxies not present in DR4). Similarly, if the small num- 
ber of colou r-based masses are neglec ted and we just use the DR4 
masses from Kauffm ann et al. (2003b) our results are unchanged. 

2 we use modelmag to measure colours throughout 

3 The scatter in the colour-estimated mass is slightly smaller using magni- 
tudes and colours transformed to z=0 rather than z=0.1 (even though the 
latter involves a smaller correction and s hould be less uncertain). T his may 
be in part due to the way in which the IKauffmann et alj fc003bh masses 
are calculated. Regardless, we present the relation in terms of z=0 values, 
which should also simplify comparison with other works. 



2.2 Line fluxes 

Line fluxes are taken from the DR7 version of the MPA database 3 
and the measurement method is described in Tre monti et alJd2004 ). 
Briefly, the spectra are compared wi th a suite of high-resolut ion 
stellar population synthesis models ( Bruzual & Chariot] 120031) to 
fit the continuum shape, thus accounting for weak metal absorp- 
tion under forbidden lines and Balmer absorption. After accounting 
for absorption and removing any remaining sky residuals, Gaus- 
sian profiles are simultaneously fitted to the emission lines, requir- 
ing that all the lines in a given class (Balmer, forbidden) have the 
same width and velocity offset. This method is therefore benefi- 
cial for measuring weak lines, since it uses constraints from multi- 
ple/stronger lines. In this paper, we will primarily be interested in 
the lines Ha A6563 (hereafter Ha) and [O II] AA 3726, 3729. We 
sum the fluxes of the two components of the [O II] 3727 doublet, 
adding their uncertainties in quadrature, and hereafter refer to this 
flux simply as [O II]. 

The MPA group derived empirical corrections to the line 
flux uncertainties via repeat observation^]. Ha error estimates are 
scaled by a factor of 2.473 and [O II] error estimates are scaled by 
a factor of 2.199. These error estimates are used when selecting 
galaxies, as we will essentially construct Ha- and [O II]-selected 
samples by requiring S/N > 2.5 in the corresponding emission 
line flux, for each of these samples. 

The SDSS is a fibre-based spectroscopic survey and thus, all 
spectral measurements are based on the flux through the fibre aper- 
ture. In order to correct these quantities to something approximat- 

4 http://www.mpa-garching.mpg.de/SDSS/DR7/ 

5 see http://www.mpa-garching.mpg.de/SDSS/DR7/raw_data.html 
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ing total fluxes, we estimate aperture corrections by comparing an 
estimate of each galaxy's total magnitude (as given by the Petrosian 
magnitude) with the corresponding magnitude measured through 
the fibre. The correction factor in magnitudes is thus given by the 
difference of these two measurements. We use the fibre and Pet- 
rosian magnitudes in the u-band as this is likely to most accurately 
trace the light from star formation. We discuss the aperture correc- 
tions in more detail in Appendix lAl 

Ratios of line fluxes (which are thus insensitive to aperture 
corrections) can be used to classify the nature of the emission aris- 
ing fr o m each galaxy, the so-ca lled BPT diagram (Baldw in et al.l 
1 1 9 8 lh . iKauffmann et alj d2003al) presented a modified version of 
the BPT diagram to classify galaxies into star-forming, AGN 
and composite types. Here we adopt the IKauffmann et alT ( l20033) 
scheme to provide a spectral type for each galaxy. For the bulk of 
the analysis we will not separate galaxies by type, but will study the 
star-formation rate density of all galaxies assuming that they each 
follow the line luminosity-star formation rate relations discussed 
below. In i]6.3.3| we discuss the impact of the different classes on 
the sample. 

2.3 GALEX photometry 

UV data is taken from the Galaxy Evolution Explorer satellite 
(GALEX). GALEX observes in two bands: the near-UV (NUV) and 
far-UV (FUV) at effective wavelengths of 227lA and 1528A re- 
spectively. A GALEX-SDSS match ed catalogue is c onstru cted in a 
manner similar to that described in Budava ri et al.l J2009h . Briefly, 
the GALEX databast0 is searched for objects within the Stripe 82 
area. To remove duplicate objects from overlapping observations, 
only 'primary' objects are retained, based on the appropriate flag 
in the PhotoObjAll table. NUV and FUV magnitudes and er- 
rors are returned as well as the coordin ates of the objects. We fol- 
low the approach o f lSalim et aljj2007l) : we use the calibrated 
magnitudes which are the flux through an elliptical aperture 
scaled to twice the Kron radius for each s ource (MAG_AUTO 
from SEXTRACTOR. lBertin & Arn outs 1996) and so are approx- 
imately total magnitudes. Object detection is run separately in the 
NUV and FUV bands. If detections are present in both bands, 
we use the elliptical aperture determined from the NUV image to 
determine the FUV flux (f uv_ncat_f lux), due to the typically 
greater sensitivity in the NUV. If a FUV detection is present, but 
not NUV, then we take the aperture parameters from the FUV im- 
age (nuv.f cat.f lux). Cross-matching is performed by associ- 
ating GALEX objects (within a 0.6 degree radius of the pointing 
centre) with the closest positional match in the SDSS data within 
4". This m ethod should produce a reasonable cro ss-matched cata- 
logue (e.g JSalim etai]|200l iBudavari et al.ll2009h . We add 0.052 
and 0.026 mag to the FUV and NUV magnitude errors re spectively 
to account for calibration uncertainty JSalim et alJ|2007T). Galactic 
extinc tion is corrected following the prescription of IWvder et all 
J2007I) . The GALEX imaging covers 90% of the 275 deg 2 SDSS 
Stripe 82 area. 



3 MEASURES OF STAR FORMATION RATE 

In this section, we present the methods to derive the default mea- 
surements of SFR using: detailed fits of spectrophotometric data 

6 see http://galex.stsci.edu/casjobs/ 



to stellar evolutionary synthesis and photoionsation models; and 
simple scaling relations between the luminosity of the Ha and 
[O II] emission lines, and the luminosity of the ultraviolet contin- 
uum. As described in the introduction, these SFR proxies have dif- 
ferent limitations and sensitivities and our goal is to compare the 
accuracy with which these measure the total SFRD of the same 
volume of the local Universe. 

3.1 04 SFRxot 

Brinchmann et al. IHooi), hereafter 04, developed a method involv- 
ing detailed modelling for deriving SFRs from SDSS spectra. Two 
different approaches are used depending on whether the galaxy is 
classified as star-forming or exhibits an AGN component, based on 
the BPT diagram. 

For star-forming gala xies, the a pproa ch is to fit the spectrum 
using a suite of iBruzual & Chariot! d2003l) models to fit the con- 
tinuum emission and also to naturally account for absorption by 
weak metal lines an d Balmer lines. Emission li nes are then fitted 
using the models of ICharlot & Longhettil toll) which model the 
emission lines of integrated spectra including a physically moti- 
vated model for dust attenuation. By fitting to grids of such models 
and calculating likelihoods using a Bayesian approach, 04 hence 
derive best fit values and associated probability density functions 
(PDFs) for various parameters for each galaxy such as SFR, metal- 
licity, total dust attenuation, etc. See 04 and references therein for 
details. This fitting technique largely depends on the luminosity of 
the Ha line, but, in practice, is estimated using all available emis- 
sion lines. 

For galaxies which are either of too low S/N to use the emis- 
sion line fitting technique, or for which the emission is from sources 
other than star-formation (i.e., AGN), an alternate approach is used. 
The relation between the D4000 spectral index and the specific star- 
formation rate (SFR/ A/*) is calculated for the star-forming sample 
just described, and this relation is then used to infer a star-formation 
rate for the remaining objects, based on a measurement of their 
D4000. 

In addition, 04 adopt an empirical, colour-dependent aperture 
correction which attempts to correct for light missed by the fibre 
using the resolved colour from SDSS photometry. Throughout this 
work, we adopt a simple aperture correction as described in ij2.2l for 
all our emission line measurements. We refer to the 04 SFR derived 
from either the emission line fitting (SFR e ) or D4000 (SFR d ) 
converted to a total SFR in this way as SFRxot ■ 

The advantage of the 04 method is that it allows detailed error 
estimates to be made by propagating full likelihood functions for 
various parameters over which the models are marginalised. The 
disadvantage is that the estimated parameters are, by their very def- 
inition, highly model-dependent and also that this technique is not 
currently widely used, particularly at high redshift where spectra 
are typically of much lower S/N. Nonetheless, we are interested in 
exploring how these complex estimates of SFR compare with more 
traditional measures. Several recent works have used the 04 SFR 
esti mates as a reference against which to compare o ther estimates 
(e.g jTrever et al.l2007llArgence & Lamareillel2009l) . and here we 
do likewise. 

3.2 HaA6563 SFR 

Taking the Ha line flux from the MPA catalogue (which has been 
corrected for absorption by the continuum fitting procedure used), 
the SFR can be estimated from the Ha luminosity using 
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SFR 

M yr- 



10 



0.4A Ho 



L(Ea) 



1.5 1.27 x 10 41 ergs- 1 



(2) 



Th e conversion fact or corresponds to the Salpeter IMF calibration 
of iKennicut multiplied by 1.5 to convert to the Kroupa 

(2001) IMF (e.g., B04). This factor is s imilar to that obtained from 
the P EGASE stellar population code dFioc & Rocca-Volme range 
1 19971) for solar metallicity, although in reality this factor depends 
on metallicity and the extinction of Lyman continuum photons be- 
fore they ionise any nebular gas. The additional coefficient assumes 
Aho mag of extinction at Ha. 

We estimate the value of extinction at Ha on a galaxy-by- 
galaxy basis using the observed B aimer decrement, the ratio of 
the flux es, /hc/Imb- W e assume an intrinsic B ai mer decremen t 
of 2.85 dOsterbrocklll989l) and the extinction law o f lSeatodft979l) , 
renormalis ed to Ry = Ayl E(B — V) = 3.1, as is commonly 
done (e.g jjansen et all200ll) . This gives 



-4; 



2.5 



k H p/k f 



■ log 



1 fua 

2.85 f H0 



(3) 



where the coefficients kup/kHa depend on the extinction law as- 
sumed, fc(A), an d kHd/kHa = 1-48 for the Seaton law (or simi- 
larly, 1.45 for fhe lCardellietal] 19891 law). The B aimer decrement- 
corrected Ha SFR (which we shall refer to as Ha corr ) is com- 
monly used as the preferred estimate of the SFR. The efficiency 
factor in converting ionizing pho tons/SFR into Ha line lum inosity 
varies weakly with metallicity dCharlot & Longhettill200ll) (which 
may be linked to a galaxy's stellar mass via the mass-metallicity 
relation), such that metal-rich stars produce lower Ha fluxes for a 
given SFR, and thus the Ha flux would overestimate the SFR for 
high metallicity galaxies. This effect is only a factor of ~1.6 over 
three decades in mass, though (04, fig. 7). 

We adopt Ha corr as our preferred estimate of SFR, but will 
later show how this compares with the more detailed modelling to 
derive the SFR. 



3.3 u-band SFR 

The ultraviolet continuum can be used to estimate the luminosity 
from young stars and hence the SFR. However, care must be taken 
to account for the contribution from older stellar populations to the 
ultraviolet light. As a simple method for accounting for the old stel- 
lar populations, we reject galaxies on the red sequence. This should 
only reject a minimal amoun t of the total SFR D in dust-reddened 
star-forming galaxies (^5%. ISalim et al' I l2007t) . We derive a cali- 
bration, using PEGASE, between the UV luminosity and the star- 
formation rate (assuming a constant SFR over lGyr), for a Kroupa 
IMF around solar metallicity (although the calibration is similar for 
sub-solar metallicities), where UV can refer to any broad band filter 
in the UV regime (e.g., u-band, GALEX NUV, FUV, etc.). The 
shape of the UV continuum (~ 1500 — 3500A) for a constant-SFR 
galaxy is approximately flat in /„ and thus the coefficient linking 
UV luminosity and SFR is relatively independent of the filter defi- 
nition. For simplicity, we adopt an average relation of 



1.4 x 10 21 WHz" 1 



(4) 



SFRyy = 1Q o.4A 

Moyr- 1 

where Auv refers to the extinction in the UV passband chosen. We 
discuss different methods of estimating Auv in £16.21 

Observed magnitudes are k-corrected to z — 0.1, again using 
KCORRECT on the ugriz photometry. To remove the contribution 
from likely non-star forming galaxies, galaxies are separated into 



red and blue populations using th e criterion u — g = — 0.033M U + 
0.75 in rest-frame magnitude^ dPrescott et alj|20091) . and only the 
blue (star forming) population is considered when measuring SFRs. 

To examine the fraction of u-band flux from various stellar 
ages within the blue galaxy population, we consider a continuous- 
SFR 10 Gyr model from PEGASE dFioc & Rocca-Volmerangel 
Taking a single model that has had a 10 Gyr constant- 
SFR at solar metallicity, the curves show the fractional contribution 
to the total flux in a given passband for bins in stellar age. For ex- 
ample, even though the SFR has been constant for 10 Gyr, nearly all 
the Lyman continuum (<91.2nm) radiation comes from <10 Myr 
stellar ages. Thus, the Ha and [O II] indicators are sensitive to star- 
formation on these timescales. The FUV (shown at ~150nm) is 
sensitive to star-formation on somewhat longer timescales (mostly 
~10 Myrs, but up to ~1 Gyr). Now, the it-band (~325nm) exhibits 
a longer tail out to much older stellar ages, since more evolved stel- 
lar populations (if present) can contaminate the u-band flux. By 
applying a colour cut in order to reject red galaxies, we have re- 
moved galaxies which will have an even higher contribution from 
stars older than ;>1 Gyr than for the continuous-SFR population. 
This leaves us with blue galaxies for which the dominant contribu- 
tion to the u-band light is from recent star-formation. For a quasi- 
continuous SFR (the blue population), about 10% of the u-band 
flux comes from stars older than lGyr for which we do not apply a 
correction. The crucial point when we come to compare the cosmic 
average of the various indicators later is that, although the emis- 
sion lines are sensitive to star-formation on much shorter timescales 
(OOMyr) than the UV estimators < 1 Gyr, both these timescales 
are short compared with the timescale over which the global star 
formation r ate of the Universe is chang ing (e-folding time ~3- 
4 Gyr, e.g.. iBaldrv & Glazebrookl 120031) . Thus, even though in- 
dividual galaxies may have longer or shorter timescales for star- 
formation than ~1 Gyr, say (and thus show differences in their 
individual SFRs estimated from Ha versus UV, for example), the 
cosmic average will show the same behaviour in all the different 
indicators, since the timescales over which these are sensitive are 
short compared with the timescale on which the global average is 
changing. 



3.4 FUV SFR 

Observed FUV— band photometry from GALEX should be a better 
proxy for the rest-frame ultraviolet continuum produced by young 
stars than the observed u-band, since it is less affected by more 
evolved stellar populations. The disadvantage of the FUV is that it 
is more susceptible to the effects of dust extinction than the redder, 
u-band luminosity. To estimate FUV luminosities, fc-corrections to 
2 = 0.1 are derived from FUV, NUV, ugriz photometry using the 
GALEX JCCORRECT routine of KCORRECT. FUV luminosity is then 
converted to SFR using Eqn.|4] 

Again, we only consider blue galaxies when measuring SFRs 
from the FUV, th is time using the r est-frame (NUV — r) < 4 
colour criterion of lSalim et al. d2007l) . According to those authors, 
95% of star-formation takes place blueward of this colour cut, and 
thus this is an observationally efficient way to select th e galaxies of 
interest whilst eliminating possible AGN contaminants. Sali m et ail 



7 The division in Prescott et al. 1200^) is given in rest-frame magnitudes, 
i.e. z=0, so although we transform luminosites to z=0.1 to minimise the 
uncertainty in the fc -correction for SFRs, we fc-correct to z=0 when sepa- 
rating red and blue galaxies 
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Figure 2. The flux contribution versus stellar age for a star-forming galaxy 
that has had a constant SFR over 10 Gyr. Model spectra were obtained from 
PEGASE with no extinction, and the different line styles (as indicated in the 
key) indicate the results for different passbands. The UV passbands 150nm 
(dashed line) and 325nm (dot-dashed line) are sensitive to star-formation on 
short timescales (ages < 1 Gyr) but the latter exhibits a tail toward longer 
timescales, meaning that more evolved stellar populations can contaminate 
SFR estimates using 325nm flux alone. See text for details. Note also that 
the Lyman continuum flux (solid blue-green line) responsible for producing 
the Ha emission line is sensitive only to ages <10 Myr. 



d2007h provide a relation to estimate the FUV extinction using 
FUV-NUV colours (their eq. 6) and we adopt this to produce an 
extinction-corrected estimated of the FUV SFR. 



3.4.1 Dust corrections to the UV SFRs 

The dust correction to the UV (u-band and FUV) SFRs is one of the 
largest uncertainties. As described above, for each passband we will 
use three different methods: a) assuming a fixed constant extinction 
of A(Ha)=l mag, b) a value scaled from the Balmer decrement 
and c) for the FUV, a correction based on the UV slope (NUV- 
FUV colour); and for the it-band, a correction which depends on 
the SFR. The first two methods require the knowledge of the rel- 
ative extinction between that of the Ha emission line and the UV 
wavelength of interest for t he stellar con t inuum . One commonly 
used prescription is that of ICalzetti et al.l {2000). This is derived 
from starburst galaxies and may not directly apply to less vigor- 
ously star-formin g galaxies dBuat et al.ll2002f) . The standard imple- 
mentation of the ICalzetti et al.l 120001) method require s the use of 
a Milky Way-like screen law (e.g. JCardelli et aill 19891) to estimate 
th e reddening suffered by the emission line. The obscuration curve 
of ICalzetti etai] d2000h is then used to estimate the reddening ex- 
perienced by the stellar continuum at the wavelength of interest 
(here, 3220A and 1390A for the z=0.1 u-band and FUV respec- 
tively), with an additional factor of 0.44 to account for the greater 
reddening experienced by the emission lines. This leads to relations 
of A(u) = 1.19 A(H a) and A(FUV) = 1.96A(Ha). Now, the 
reason for using the two different obscuration curves (foreground 
screen for lines and an embedde d geome t ry for the stars) for a 
starburst galaxy is sketched out in ICalzettil ( BOOlh : but is likely an 
overestimate for more modestly star-for ming galaxies. In or der to 
attempt to mitigate this effect, we use the ICalzetti etai] d2000h ob- 
scuration curve for both the stars and the emission lines (but again 



use the factor of 0.44 to convert between the relative reddening in 
the two). This yields relations of 



A(u) = 0.87A(Ha) 
and 



(5) 



A(FUV) = 1.43A(Ha). (6) 

This a pproach is somewhat similar to the model of ICharlot & Falll 
(2000) who use a single attenuation law for birth clouds and the 
diffuse ISM, with a different normalisation between the two. These 
numb ers are also broad ly in agreement with the obs ervational re- 
sult oflBuat et al.l d2002l) who found that the standard Calze tti et al.l 
( 2000) prescription overestimates the UV extinction by ~0.6 mag 
when applied to non-star bursting galaxies. We obtain broadly sim- 
ilar results (A(u) = 0MA(H a ) and A( FUV) = 1.55A(H a ))if, 
instead of using the Calzet tiet al.l J2000h obscuration curve to cor- 
rect both the stars and the gas, we use only the Milky Way extinc- 
tion curve of lCardelli et alj fl989l) for both (again using the factor 
of 0.44 to convert between stars and gas). 

Regardless of the exact motivation, we will show that these 
scalings between extinction at Ha and the UV (eqn. [5] & eqn|6} 
yield reasonably good agreement between the SFRD estimated in 
the UV and Ha and so we adopt these as our default normalisa- 
tions. Our SFRDs may straightforwardly be renormalised to as- 
sume a different dust prescription by simply multiplying the SFRD 
by iQ 0A (° ld - new ) where old represents the coefficient in eqn.[5]or 
[6] and new is the desired alternate coefficent. We will discuss the 
choice of extinction correction further in ! 



3.5 [OII]A3727SFR 

The SFR measured from [O II] luminosity may be estimated by 
scaling between the [O II] and Ha luminosity using 
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(7) 



where rii ncs is the ratio of extinguished [O II] to Ha flux. In 
the absence of bette r information, a ratio of « 0.5 is typi- 
cally assumed (e.g.. lKennicutt| [T998). We take this value of 0.5 
and assume 1 mag of extinction at Ha (as is commonly done) 
as our nominal [O II] SFR. This gives a nominal SFR relation, 
SFRo/iMQyr' 1 ) = L([OII])/(3.80 x 10 40 ergs' 1 ). 

It is well known that the r atio, rij ne s, depends on metallicity 
and ionisation parameter (e.g.. iKewlev et al" I l2004l and references 
therein). As a default model, we assume a fixed conversion between 
[O II] luminosity and SFR, as is often done. In ^4]we will consider 
empirical corrections which attempt to account for this (i.e., with- 
out needing to assume a fixed rii nes ). 



4 EMPIRICAL CORRECTION TO THE [O II] SFR 

The [O II] luminosity from the [O II]AA3726, 3729 doublet is less- 
directly coupled to the ionising flux from young stars than Ha, and 
thus suffers larger uncertainties as a SFR indicator. As mentioned 
above, many of the systematic uncertainties in the SFR indicators 
may be related to stellar mass (e.g., through the mass-metallicity 
relation). So, it is of interest to examine how the relation between 
[O II] SFR (estimated simply from [O II] luminosity with a con- 
stant conversion to SFR) compares with a more direct SFR indi- 
cator, on a galaxy-by-galaxy basis, as a function of stellar mass. 
Fig. [3] shows [Oil] SFR versus the Balmer decrement-corrected 
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Figure 3. Comparison of Balmer decrement-corrected Ha SFR against constant extinction [O II] SFR rate separated into bins of stellar mass. Note: fit is 
performed over full range of SFRs, but plots are scaled to show most populated regions for clarity. Best fit lines with lcr uncertainties are shown in green. 
Dashed blue line indicates the slope adopted in the best-fit fitting function in each mass bin. The solid black line indicates equality, for reference. The final 
panel shows the best fit ratio of S F R,\q j j^l S F Ru a as a function of stellar mass (black points with error bars); the green line is the adopted tanh correction. 
For reference, blue diamonds indicate the ratio of the total SFRD from [O II] to SFRi* t (Fig. [§}, but this is not used in the fit. See text for details. 



Ha SFR. For this comparison, only galaxies identified as belong- 
ing to the star-forming class by the BPT classification described in 
$2.2\ are considered. In each mass bin, a robust linear regression 
fit (required to pass through the origin) is performed rejecting out- 
liers. The best fit and its la uncertainty are indicated. It is clear 
from the figure that the slope of the relation is a strong function 
of mass. In order to model this dependence, the best-fit slope and 
its uncertainty in each mass bin is fitted with a tanh function of 
mass (shown in the final panel of Fig. [3}. The tanh function effec- 
tively captures the shape of this relation, being approximately hor- 
izontal at high and low masses and switching between these two 
levels in an approximately linear fashion at intermediate masses. 
The empirically-corrected SFR, SFR ernp , corr is obtained from the 
nominal [O ITJ-SFR (EqnEJl, SFR , 

WB _ SFRp 

J r riemp cor , (o) 

a tanh[(x- — b)/c\ + a 

where x = log(M»/M ) and a = -1.424, b = 9.827, c = 0.572, 
and d = 1.700, and SFR / {M @ yr" x ) = L([OII])/(2.53 x 
10 40 erg s" 1 ) from above. This, should give the most reliable 
estimate of SFR from [O II] luminosity. The correction between 
[O II] and Ha corr SFR encodes all the information about the 
metallicity and dust dependence of [O II] luminosity on SFR in 
each mass bin. In principle, this can now be used to correct 
[O II] luminosity to SFR in higher redshift. If this correction is used 
at higher redshift, this implictly assumes that the mass-dependence 
of extinction, metallicity, etc. (which are all encoded in the correc- 
tion) behave in the same way at higher redshift as locally. We note 
that this (possibly incorrect) assumption is commonly used anyway 



at higher redshift when locally-derived relations (e.g., Eqn. O are 
applied to higher redshift galaxies. We will consider this point fur- 
ther in § 16.41 but for the moment note that the factor most likely to 
evolve between here and the higher redshift Universe is likely to be 
the amount of dust in a galaxy of given mass. Naively, one might 
assume that the dust content of a galaxy would increase from high 
to low redshift, as successive generations of stars pollute their host 
galaxy. However, it is equally possible to make simple arguments 
for the contrary evolution based on dust geometry, for example. To 
simplify comparison with future work, Fig.|4]shows the amount of 
dust extinction, as measured from the Balmer decrement, as a func- 
tion of stellar mass. The solid line shows a fit to the mass-dependent 
extinction, which may be modelled as 

A Ha = a + 61og(M*/M ) + c[log(M,/M )] 2 (9) 

with the values a = 51.201, b = -11.199, c = 0.615 and setting 
the fit to a constant value below log(M»/M )^9.O. For reference, 
the dashed line shows the average mass dependence of extinction 
found by 04, taking the value from the peak of their likelihood 
functions in each mass bin from their fig. 6. The coefficients are 
approximately a = —3.39, b = 0.46 (c = 0). Their relation is dif- 
ferent from ours since we use the Balmer decrement whereas they 
estimate the dust by simul taneous fits to many parameters in the 
ICharlot & Longhettil J200lh models. These models include a com- 
ponent of diffuse nebular emission which has a lower effective tem- 
perature than the canonical value assumed for HII regions and, in 
effect, leads to the adoption of a non-constant ratio for the Balmer 
decrement (J. Brinchmann, priv. comm.). This is phrased as "metal- 
licity dependence of the Case B Ha/H/3 ratio" in 04. It can be seen 
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Figure 4. Extinction in the Ha line (estimated from the Balmer decre- 
ment) as a function of stellar mass. Points show measurements for each 
galaxy, and the solid line is our quadratic fit described in the text. Dashed 
line shows the relation inf erred from 04 (which uses fits to the models of 
ICharlot & Lo nghetti 2001 to attempt to estimate the dust). Larger filled cir- 
cles show extinction scaled from IR estimates, discussed in i]6.2.4l 

from Fig.[4]that, for low mass (log(M, /Mq)~9) galaxies, allow- 
ing for a non-constant Balmer decrement ratio leads to extinction 
estimates approximately twice as large as those inferred from our 
constant Balmer ratio estimate. When estimating SFR, this is al- 
most precisely cancelled by a corresponding v ariation in conversion 
factor between Ha luminosity and SFR in the lCharlot & Longhettil 
models (fig. 7, 04). Given the uncertainties in fitting such 
a diffuse component, we prefer to adopt the more traditional con- 
stant Balmer decrement and constant conversion from Ha lumi- 
nosity to SFR. 04 showed that the these two approaches give very 
similar estimates of SFR (and we will also show that this is true 
when we compare SFRDs derived from the different indicators). 
Furthermore, we also show in Fig. [4] (larger filled circles) an al- 
ternate estimate of extinction based on FIR measures (discussed in 
S]6.2.4t which can be seen to be in better agreement with that de- 
rived from our constant Balmer decrement estimate than the vari- 
able ratio model fits. 

These various parameterisations of how dust extinction varies 
with mass will allow the reader to remove the dust term from our 
empirical correction and substitute an alternate correction, if de- 
sired. In addition, Fig. [4] in conjunction with Fig. [3] shows how 
much of the systematic offset of [O II] to Ha SFR is due to the 
effect of dust. 



5 EMISSION LINE LUMINOSITY FUNCTIONS 

Before proceeding to more detailed studies using SFRs based on 
our [O II] and Ha measurements, it may be instructive to compare 
the distributions of our measurements with those of other surveys. 
The simplest useful statistic employed by many works is that of the 
luminosity function. In order to more easily compare with other 
works we use our observed quantities directly. These may simply 
be rescaled using the relations given above, if the reader desires to 
see these in terms of SFR functions, rather than observed emission 
line luminosity. Another useful result of such a comparison is that 
it will allow us to confirm the depth of our emission line flux limit, 

by examining the faintest line s visible. 

We follow the method of lzhu ct al. (2008) and use a 1/1 ma . r 



Figure 5. Luminosity functions for Ha and [O II] in Stripe 82. Red and 
blue lines with error bars show observed (i.e., no correction for extinc- 
tion) luminosity functions for Ha a nd [O II] respectivel y. For compari- 
son, observed [O II] (z<0.045) from Gall ego et all 120021) (blu e squares); 
extinction corrected Ha (z<0.()45) from Galleg o et al] j 19951) (red aster- 
isk s); non-extinction corrected Ha from the z < 0.01 narrowband survey 
of (James et al. 2008) (P. James, priv. comm.); and observed Ha from the 
z~0.08 narrowband survey of Ly et al. (2007) (red triangles) are shown. 
Double power law fits to the Stripe 82 data are overplotted. Dotted ver- 
tical lines (red: Ha; blue: [O II]) illustrate luminosities corresponding to 
SFRs?» 0.1M©yr for the simple constant conversions. See text for dis- 
cussion. 

technique to estimate the number of objects in bins of log luminos- 
ity (see also ^6. H . Error bars are estimated u sing scaled Poissonian 
distributions, as described in Izhu et al.l ( 120081) . The resulting LFs 
are listed in Table Q] and plotted in Fig. [5] Filled red circles show 
the Ha LF and open blue diamonds show that for [O II]. Neither 
LF is well fitted by a Schechter function, which is the normally- 
employed fitting function for galaxy LFs. As noted by Izhu et al.l 
(2008) at z~l, a better function to use is a double power law. 
They found that the higher redshift [O II] LFs are well-described 
by a power law with a break around log Lou ~ 42 with a slope, 
a ~ —3 (dN/dL oc L a ) at higher luminosities, and a ~ —1.3 
faintward of the break. Fitting a similar function, we find that the 
Ha LF is well-described by a function with ctbright — —2.48, 
otfaint = —1-27 and a turnover luminosity, log Lto = 41.33. For 
the [O II] LF, these values are airight = —2.83, a faint = —1.67 
and log Lto = 41.03. The dotted vertical lines Ha and [O II] (red 
and blue respectively) show the median luminosities at which the 
S/N falls to 2.5 (our selection limit). 

Local LFs from the literature are overplotted on the figure. 
There are few recent [O II] LFs, since Ha is the preferred SFR in- 
dicator locally, and is more easily accessible (at very low redshifts, 
few spectrographs have sufficient blue sensitivity to the observed 
wavelength of [O II]) . The observ e d [O I I] (i.e., uncorrected for 
extinction) LF from Galle go et al.l d2002l) (converted to our cos- 
mology) for a sample of ~ 150 galaxies at z<c0.045 is shown as 
open blue squares. This is in excellent agreement with the Stripe 
82 measurement at the bright end. The turnover occurs in the 
same place in both datasets. Howev er, the faint en d slope appears 
somewhat shallower than ours in the lGallego etall d2002l) sample. 
This may be a sign of incompleteness in the latter dataset, and 
the much greater size and un iformity of the SDSS data should be 
favoured. iGallego et af] d 19951) presented an Ha LF from the same 
survey as their [O II] LF. However, in this case they only pub- 
lished extinction-corrected Ha luminosities, shown as red asterisks 
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Table 1. Luminosity functions for [O II] and Ha in Stripe 82, plotted in 
Fig. \5\ Number densities, <fi, are observed values with no corrections for 
incompleteness or internal extinction. Incompleteness will begin to set in 
fainter than log L R240.3. See text for details. 
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in Fig. [5] Even so, these are in excellent agreement with our ob- 
served (non-extinction corrected) LF (except for their faint est bin, 
which is likely again incomplete). At lower luminosities. iLv et al.l 
d2007h published an Ha LF based on narrowband imaging of the 
Subaru Deep Field at z~0.08. This was drawn from too small an 
area to constrain the bright end of the LF, however the faint end 
appears to have approximately the same slope (a ~ —1.6) as ours, 
but with higher normalisation. This may be due to cosmic vari- 
ance, due to the small field size, or possibly some other system- 
atic difference between the datasets based on narrowband imag- 
ing and flu x-calibrated spectr oscopy. The narrowband Ha imaging 
survey of I James et al.l d2008l) similarly shows a steeper faint end 
slope (fainter than our completeness limit) suggesting that this is 
indeed likely due to completeness effects. Brighter than this, the 
Ijames et"al] d2008l) Ha luminosity function is in good agreement 
with our measurements for all but the brightest point which dis- 
agrees at ~2er. This could again be due to cosmic variance, due to 
the relatively small area, or to evolution, since this sample is se- 
lected to be z < 0.01. 

All our results for Ha and [O II] so far have considered galax- 
ies with S/N>2.5 in the given emission line, since we require a 
minimum S/N in each emission line in order to assign an object an 
SFR from the emission line measurements. An important question 
for understanding the SFRD measurement is how well an [O II]- 
selected survey (or Ha-selected survey) mimics an SFR-selected 
survey. The above comparison of the [O II] and Ha LFs in Stripe 82 
shows that the number densities of Ha-selected and [O II]-selected 
objects at a given luminosity (SFR), close to the limiting luminos- 
ity, are comparable. In order to more effectively study the effect of 
selection, we next consider directly the star formation rate density 
(SFRD) and its mass dependence. 



6 THE STAR-FORMATION RATE DENSITY 

Whilst the comparison of SFR indicators on a galaxy-by-galaxy 
basis may be the most instructive method for understanding sys- 
tematic differences between indicators, perhaps a more useful way 
for understanding results from large galaxy surveys is to exam- 



ine the ensemble properties of galaxy populations such as the star 
formation rate density (SFRD). The advantage of studying this 
cosmically-averaged quantity is that the timescale sensitivities of 
the different indicators considered are all short compared with the 
timescale on which the cosmic average is changing, and this allows 
a fair comparison between the different indicators. Indeed, a com- 
parison of UV SFR versus Ha SFR for individual galaxies might 
be expected to show significant scatter, since these indicators are 
sensitive to star-formation on different timescalejf], but averaged 
over large galaxy populations, such differences are reduced. 

6.1 Calculating SFRD 

The star formation rate density in the m th m ass bin, psfr d, is 
calculated using the 1/V m ax methodology (e.g.. lFeltenll 19761) 

1 \ - SFRj 

PSFR, m = -r-j T7 } J- (10) 

A log M ^ Vmax.i &i 

where SFRi is the SFR of the i th galaxy, V ma x,i is the maximum 
volume out to which it could be observed in the survey, and the 
weighting, Cj, is the spectroscopic completeness described in ^2] 
If a galaxy is bright enough to be seen over the full redshift range 
(0.032 < z < 0.20) down to the magnitude limit of the survey 
(r = 19.5), then it would be drawn from a volume of 1.5 x 10 7 
Mpc 3 . For galaxies which could not be seen out to the high redshift 
limit, it is necessary to calculate the redshift at which it would be- 
come too faint to be included. This procedure involves knowing the 
spectral type of the galaxy in order to estimate its fc-correction (at 
its maximum redshift, again using KCORRECT), and ultimately the 
maximum redshift (and hence Vmax) must be calculated iteratively. 

Fig. |6] shows the SFRD as a function of stellar mass for the 
04 SFR Tot . For reference, the 04 estimate of the total SFRD from 
a larger area of SDSS (DR2, and in the redshift range 0.005 < 
z < 0.22) is overplotted as the green dashed line. Poisson uncer- 
tainties based on the number of galaxies in each mass bin are much 
smaller than the systematic uncertainties and so error bars for each 
individual method are omitted from the plot for clarity. 

6.2 Comparison of indicators 

6.2.1 Emission line indicators 

Fig. |6] shows the comparison between the Ha SFRD (left panel) 
measured using the Balmer decrement-correction, Ha corr (red 
solid line), the average mass-dependent extinction using our re- 
lation from Fig. [4] (red dashed line) and a constant correction of 
Ah a = 1 mag (red dotted line). It is clear that the various different 
prescriptions (including simpler aperture corrections and using the 
same indicator for all galaxies regardless of spectral class) for Ha- 
inferred SFRs produce reasonable agreement with the modelling of 
04, particularly Ha corr . The biggest discrepancy is for the constant 
extinction method, as 1 mag of extinction at Ha is clearly an over- 
estimate for low mass (log(M»/M0)~9) galaxies, and an under- 
estimate for higher mass (log(M*/MQ)>10), as can also be seen 
from Fig. [4] The mass-dependent extinction correction is an im- 
provement over the constant correction model, but overestimates 
the SFRD at masses above the peak. We note that if we impose 

8 We only use the galaxy-by-galaxy comparison to correct the [O II] SFR 
to that of Ha since both these indicators are essentially sampling the same 
stars. 
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Figure 6. The SFRD as a function of stellar mass for the different emission line measurements of SFR. The left panel shows rates based on Ha (red lines) : solid 
line is Balmer decrement-corrected (on a galaxy-by-galaxy basis) Ha; dotted line is for 1 mag of constant extinction at Ha (in Eqn.[2); dashed line assumes 
average mass-dependent extinction from Fig. [4] The solid black line shows the SFRD calculated using SFR^ot (the modelling technique of 04); and the 
dashed green histogram shows the SFRD from 04 for the larger SDSS DR2 sample. The right panel shows the SFRD from [O II] assuming 1 mag of constant 
extinction at Ha (dotted blue line) and that using the empirical correction to [O II] (solid line). The empirically-corrected [O II] SFRD only underestimates 
the Ha SFRD due to the [O II] flux limit (see also Fig.[9|and text for discussion). For reference, Ha CO rr and SFR^ ot results are repeated from left panel. 
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Figure 7. The total SFRD as a function of limiting SFR in bins of stellar mass for SFR^ ot (black solid line), Ha corr (dashed red line), and [O II] emp , corr 
(dot-dashed blue line). Units of SFRD are 10 — 3 Mq yr -1 Mpc~ 3 . In stripe 82, the SFRDs based on SFR^ot and Ha corr largely appear to have converged 
for a limiting SFR of log(SFR) ~ — 1. However, the [O II] based SFRDs appear to not have converged in the highest mass bins, even below this limit. See 
text for discussion. 



an upper limit of Ana = 2.0 mag on the mass-dependent extinc- 
tion, this overestimate of the SFRD at higher stellar masses disap- 
pears, and thus it may just be that the fit in Fig. [4] is asymptoting 
to too high a value of the extinction. The dashed green histogram 
shows the SFRD of 04, measured from the larger area DR2. This 
is very slightly different from the result measured in the Stripe 82 



area using their SFR*Tot indicator (black solid line). The peak of 
the former is offset to slightly lower masses (by ~0.2 dex). This 
may be due to evolution (see next section); the 04 sample extend- 
ing to slightly lower redshift; cosmic variance; or possibly differ- 
ences in constructing the SFRD histogram. 04 use the full proba- 
bility distribution functions for the SFRs of each galaxy measured, 
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Figure 8. As Fig. [6] but for the photometric SFR indicators. Blue lines in the left panel show u-band derived SFRDs. Solid line is corrected for extinction using 
the Balmer decrement, dotted line assumes constant extinction (equivalent to 1 mag at Ha), and dashed line assumes SFR-dependent extinction. Right panel, 
blue lines show measurements from FUV. Dashed line indicates extinction-corrected SFRD using the UV slope (NUV-FUV colour), dotted line assumes 
constant exti nction (1 mag at H a), and dashed line assumes SFR-dependent extinction. Dashed histogram shows the (approximately FUV)-derived SFRD 
measured by ISalim et alj2007i) by SED fitting. 



whereas we have summed discrete values to construct the his- 
togram. Regardless of the cause, the difference between the Stripe 
82 SFRTot SFRD and the 04 result is small compared with the 
systematic offsets between indicators, which are of interest here. 

The right panel of Fig.[6]shows the comparison for the [O II]- 
based SFRDs (blue lines). It is clear that the constant-extinction 
[O IIJ-SFR (blue dotted line) is a very poor tracer of the SFRD. 
The empirically-corrected [O II] SFRD traces the Ha corr and 
SFR To t SFRDs for masses up to log(M,/M )~ 9.5. Thereafter 
it begins to systematically underestimate the Ha-based SFRDs, 
but still peaks at almost the same mass. (In S |6,3.1| we will show 
that this discrepancy is simply due to the [O II] flux limit, and the 
agreement between the Ha corr SFRD and that of the empirically- 
corrected [O II] SFRD is excellent.) In order to quantify the agree- 
ment between different indicators, we consider three measurements 
from the SFRDs: the position of the peak, the integrated SFRD 
within 8.5 log (M./M©)^ 12.0 and the integrated SFRD within 
8.5 ^log(Af*/M0)^ 10.4 (since many indicators begin to di- 
verge above this approximate position of the peak). These values 
are listed in Table [2] 

Before discussing further the differences between the SFRDs 
using different estimators, it is pertinent to consider how the lim- 
iting SFR probed by each indicator affects the final SFRD in each 
mass bin. i.e., whether the SFRD has converged. Fig. [7] shows the 
cumulative SFRD as a function for SFRTot (solid black line) and 
the best estimators for each of Ha (dashed red line) and [O II] (dot- 
dashed blue line). The emission line estimators are plotted down to 
the limiting SFR, where this is defined as the median SFR which 
gives a S/N of 2.5 in the relevant emission line. In addition, the 
FUV SFR corrected using the UV slope is included (green dotted 
line), since it does not rely on any emission line estimator (nor on 
their associated aperture corrections, indeed not on any spectro- 
scopic measurement except redshift) and is thus independent of the 
other methods. 

For the lowest stellar masses, all methods except the FUV give 
very similar answers. As well as reaching comparable total SFRDs, 
each of these methods turn over at comparable SFRs (log(SFR) ~ 
-1). Above a mass of log(M*/M Q )~10, the estimated SFRDs 
all begin to converge to similar totals, however, towards increasing 



masses, it can be seen that the [O II] estimate is still rising as a 
function of limiting SFR, i.e., the [O II] SFRD does not converge. 
This can be understood in terms of Fig.[3](most clearly from the last 
panel). At these stellar masses, the ratio between SFR (as traced 
by Balmer decrement-corrected Ha) to [O II] constant luminosity 
SFR is ~ 10: 1 . This means that [Oil] is approximately an order 
of magnitude less sensitive to SFR than Ha and becomes less so 
as stellar mass increases. Thus, to the limit of the SDSS Stripe 82 
data, the [O II]-inferred SFRD has not converged to the Ha value, 
for the highest mass galaxies. For Ha, SFRTot, and FUV, all es- 
timators show turnovers in the SFRD at similar values of SFR, but 
the final totals achieved are different. This implies that all values 
have converged, but to different limits reflecting systematic differ- 
ences in the SFR estimators. These will be discussed in the context 
of the mass-dependent SFRD. 

6.2.2 UV indicators 

Fig. 00 shows the same plots as Fig. [6] for the photometric (u-band 
left panel and FUV right panel) SFR indicators. To the methods 
already described in |J3] to measure photometric SFRs, we add 
the following methods to attempt to correct for dust extinction. 
A u , is approxim ated using the mass -dependent extinction relation 
(Eqn.[9j and a lCalzetti etail feOOCl) extinction law; and also cal- 
culated using the Balmer decrement to calculate the extinction at 
Ha, as described ab ove, and transforming this to 3270A using the 
ICalzettietal1 f2000) law. For the small subsample of galaxies with 
insufficient S/N to measure the Balmer decrement, we estimate the 
extinction at Ha from Eqn. [9] As a further possible improvement 
to the u-band SFR, we consider the SFR-dependent e xtinction cor- 
rection. We use the equations of Hopk ins et al. I fcOOlh to determine 
a relation between extinction and the observed UV luminosity. The 
SFR was then determined using Eqn. [4] In addition, we use the 
Balmer decrement as described above to estimate the extinction at 
1400A and provide an alternate dust-corrected estimate of the FUV 
SFR. 

The biggest discrepancy between indicators for low stellar 
masses occurs for the FUV SFRD compared with the other meth- 
ods. From Fig. [8] it can be seen that common to both sets of UV 
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data, the estimators using only photometry (and photometric es- 
timates for the dust correction) show an excess SFRD over the 
emission line based estimates at low stellar masses. The UV slope- 
corrected FUV SFR (i.e. that plotted in Fig. Q uses the approxi- 
mate calibration of dust extinction from the UV colou r calibrated 
from the more detailed modelling of lSalim et al. U2007h . Those au- 
thors showed that this approximation has a tendency to overesti- 
mate the SFR (compared with their SED modelling) at low values 
of SFR, and underestimate it at high val ues. However, the dashed 
histogram shows the SFRD derived by ISalim et al. I d2007l) using 
their modelling, which is in very good agreement (certainly at the 
low mass end) with our estimate. Indeed, a galaxy-by-galaxy com- 
parison with Ha corr (similar to Fig. [3) shows excellent agreement 
between the FUV slope-corrected SFR and Ha corr with a small 
systematic underestimate of the latter by the former of ~10% but 
little mass dependence until \og(M, /Mq)> 10.5, when the ratio 
drops further. This is better than might be expected for a galaxy- 
by-galaxy comparison, since the different timescales to which the 
UV and Ha are sensitive would be expected to introduce signif- 
icant scatter. Thus the result that the low mass end of the SFRD 
is predicted to be higher relative to the other methods when using 
UV data appears to be robust. The M-band and FUV SFRDs using 
the B aimer decrement to estimate the dust are in closer agreement 
with the other estimators (which also largely or entirely use the 
Balmer decrement to e stimate extinction) . Since our FUV SFRD 
agrees well with that of Sal im et aT] ^20071) and our Hq eorr i s close 
to that of 04, we refer the interested reader to Sali m et al. I J2007l) 
for a discussion of the likely reasons for the discrepancies between 
these estimators. They also found that the main disagreement is due 
to the way dust is estimated in the two methods. They went through 
an extensive analysis trying to reconcile the dust estimators, and 
showed that the SED-fitting (i.e., UV-slope based) extinction esti- 
mate is noisier than that from the Balmer decrement, and thus likely 
overestimates the extinction at low stellar masses. 



6.2.3 FIR indicators from the literature 



cresting to co i 

tion derived from Salim et al. U2007l) 's fitting to UV+optical SEDs 
with a similar relat ion derived from a different SFR indicator. 
ICortese et al.l d2006h estimated UV extinction by comparing UV 
with total IR luminosities for a local sample of late-type galax- 
ies with GALEX (UV) and IRAS (FIR) photometry. The combi- 
nation of UV and FIR should sample both the unobscured and 
obscured components of star-formation , and so, when appropri- 
ately combined (e.g., Martin et al. 20071 eqn. 7), should give the 
best estimate of total SFR for a galaxy. This ratio is usually de- 
noted the Infrared excess (IRX = \og{Lm/ Luv)- Taking the re- 
lation betwe en IRX and the ex tinction in the GALEX FUV-band, 
Af uv , fromlBuat et alj 1 20051) and combining this with eqn. 4 & 
5 from lCortese et alj feOOfj) . we obtain a relation between Afuv 
and (FUV-NUV) colour. This relation, based on IR estima tes of the 
extinction, compares very well with the lSalim et alj J2007I) relation, 
with an rms difference of 0.12 mag over the range of (FUV-NUV) 
colours of star-forming galaxies in our sample. We can approximate 
the relation as A FUV ~ 2.61(fW - N UV) + 0.409 which 
may be compared with I Salim et alj J2007l) 's eqn. 6 of Afuv = 
2.99(FUV - NUV) + 0.27 for galaxies passing their blue colour 
cut. Thus, our FUV estimate of the SFRD using the UV slope to 
estimate the extinction is relatively robust to the method by which 
the extinction is derived. 

Therefore, using the Balmer decrement to correct the FUV 



Table 2. Properties of the SFRD for the various different SFR in- 
dicators. Columns are: method used; log(M*/MQ) at which the 
SFRD peaks; integrated SFRD over the widest mass range considered 
(8.5 <log(M*/M )< 12.0); integrated SFRD from the complete- 
ness limit to the approximate peak (8.5 <log(Af*/Af0)< 10.4) - 
this latter range is chosen since the largest disagreement occurs at mass 
higher than the peak SFRD. Integrated SFRDs are given in units of 
X 10 — 3 MQyr~ 1 Mpc~ 3 . (* - histogram is bimodal, this value is midway 
between the two peaks. ** - exact normalisation depends on dust conversion 
from Ha to UV, see gXT) 



Method 



peak mass 



fs> dM 



Js.5 P dM 



SFR Tot 


10.6 


15.24 


5.86 


Ha Balmer deer. 


10.6 


13.41 


5.65 


Ha(const.) 


10.0 


12.94 


7.79 


Ha mass dep. 


10.4 


16.11 


7.92 


[O II]emp,corr 


10.4 


10.58 


5.18 


[O U]MB,corr 


10.0 


7.19 


4.90 


[O II](const.) 


9.8 


13.25 


10.61 


u Balmer deer.** 


10.8 


15.16 


7.06 


u mass dep.** 


10.4 


18.53 


9.49 


u SFR-dep. 


10.4 


15.45 


8.23 


u (const.)** 


10.0 


13.36 


8.15 


FUV UV slope 


10.4 


16.21 


8.78 


FUV Balmer deer.** 


10.4* 


13.00 


4.24 


FUV const.** 


10.4 


9.49 


6.10 


04 


10.75 


14.78 


6.53 



SFR possibly represents the best combination of methods. The 
FUV is a good tracer of star-formation (once corrected for extinc- 
tion) and is much less affected by the presence of an AGN or the 
lack of a measurable Ha emission line, and the Balmer decrement 
likely offers the best estimate of extinction. We reiterate that for 
galaxies for which the lines necessary to measure the Balmer decre- 
ment possess insufficient S/N, we use the average mass-dependent 
extinction, and thus are not limited to requiring Ha such as in the 
Ha-selected sample. The fact that this hybrid method follows quite 
closely Ha corr reinforces the suitability of this method. 

A further point worth noting from Fig. [8] is that the u-band 
SFRD, even uncorrected for dust, provides a surprisingly accurate 
estimate of the SFRD compared with the more observationally ex- 
pensive estimators. A mass-dependent extinction correction such 
as Fig.|4](under the assumption that the locally derived value is ap- 
plicable) converted to the appropriate wavelength would provide 
an even better estimate, and thus these are useful indicators for 
higher redshift surveys. A better estimate still can be made using 
bl uer data (FUV) plu s a correction for the UV slope as described 
in ISalimet all (2007). The greater sensitivity of the bluer data to 
the dust correction can clearly be seen by comparing the constant 
luminosity estimators (dotted lines) in both panels of Fig. [8] 



6.2.4 Other empirical corrections 

iMoustakas et alj d2006l) propose an empirical correction from 
[O II] luminosity to SFR using the rest-frame B-band luminosity, 
Lb- They give a method for computing rest frame B-band lumi- 
nosities that can be calculated from SDSS g and r photometry (their 
eqn. 1 and KCORRECT) and plot the ratio between Ha-derived SFR 
and I/[on] as a function of Lb (their fig. 19). We construct Lb- 
corrected [O II] SFRs for the Stripe 82 data following this approach 
and interpolating their fig. 19 to derive a correction for each galaxy 
on the basis of its Lb - This produces the dotted curve in Fig. [6] It 
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can be seen that this correction severely underestimates the SFRD 
at high stellar masses. 

Iweiner et alj d2007h derive a correction for [O II] SFR based 
on rest-frame //-band magnitude, from observed A'-band data at 
z~l, by comparing the observed [O II] luminosity with SFR de- 
rived from Ha emission. They give this correction as 0.23 dex 
mag -1 of rest-frame //-band. We synthesise rest-frame //-band 
magnitudes from the stellar mass es and rest - frame (g — r) colours 
described in §2. II using table 7 o f lBelletal.ll2004h . We set the nor- 
malisation of the //-corrected [O II] SFR by requiring the SFRD 
at log(A/*/A/0)~lO to agree with that of [O II] C orr- This gives a 
curve (not plotted for clarity) which match es very closely over the 
whole mass range considered to that of the Mou stakas et alj J2006l) 
correction plotted in Fig.(6] and so these two indicators may be con- 
sidered equivalent. Neither of these corrections adequately capture 
the large corre ctions necessary to [O II] SFR at high stellar masses. 

Recently iPannella et alj ( 120090 estimated the median SFR in 
bins of stellar mass at z~2 using photometric redshifts to select 
galaxies and determine stellar mass, and stacked 1.4GHz radio 
emission to determine total SFR in these mass bins. They found 
that the SFR is almost a linear function of log(M»/M0), i.e. a 
power law in M* with an exponent of order unity0 By comparing 
the SFR estimated from the rest-frame 1500A luminosity with the 
stacked radio SFR (assuming the latter is an accurate measure of 
the total SFR), they derived an empirical relation for the extinction 
as a function of stellar mass (assuming the difference between the 
two SFRs is entirely due to extinction in the 1500A flux estimates): 

A 150 o = 4.O71og(M,/M ) - 39.32. (11) 

Applying this (rather steep) correction as a function of stellar mass 
(converted using our assumed extinction described in £13.4. 1 b to our 
Stripe 82 u-band estimated SFR results in a vast overproduction of 
the z~0.1 SFRD. The attenuation predicted at log(A/*/M Q )~l 1 
is over 6 magnitudes, or a factor of 250, and thus the peak of the u- 
band SFRD would be off the top of the plot in Fig.[8] Similarly, us- 
ing our extinction law to convert the 1500A continuum attenuation 
to attenuation in the gas at the wavelength of [O II] or Ha would 
lead to larg e overestimates in the SFRD for these other indicators. 
Clearly, the lPannella et alj J2009h prescription is not a good esti- 
mate for the mass-dependent extinction correction in the local Uni- 
verse. This suggests that either the dust properties as a function 
of mass have evolved strongly between z~2 and z~0, such that 
extinction is much lower now; or possibly that observational sys- 
tematic errors in the z~2 data lead to an overestimate of the correc- 
tion. If the stacked galaxy samples contain contamination in the ra- 
dio data (such as might be produced by low luminosity AGN) then 
the radio-S FR (and hence extinc tion correction) would be overesti- 
mated (seelP^eUaetaDHQOi). 

We note that other groups jMartin et ail 120071 : iBuat et alj 
1 20071 : [iglesias-Paramo et al. 200% have used the infrared excess, 
IRX, as a measure of extinction and found more modest UV ex- 
tinctions than predicted by eqn.[TT] These works do not explicitly 
give a com parable relation, fitt ed as a function of stellar mass, but 
we can use lMartin et alj J2007l) 's IRX versus stellar mass data and 
convert IRX to A(Ha) as follows. We first take t he rel ation be- 
tween IRX and UV slope, /3 from ICortese et alj yOOg, eqn. 5). 
Now, we just need to convert between UV slope, /? and A (Ha). 



9 |PunneetalJj2009l) came to similar conclusions using a similar technique 
on a different field, but with a somewhat shallower relation. 



ICortese et alj fcOOfJ) present two relations between these quanti- 
ties, one fitted from their o wn data (their eqn. 9) and the other taken 
from lCalzetti et all d2000f) (the former's eqn. 8). However, we may 
fit our own relation between j3 (simply related to the (FUV-NUV) 
colour, e.g. JCortese et alj2006l . eqn. 4) and A(Ha) taken from the 
Balmer decrement. Doing this we find /3 ~ 1.26A(H a ) — 1.79 
but the scatter is large and the fit is not well-constrained. We show 
the results of taking th e z=0.1 stellar mass -IRX data (adjusted to 
our IMF) from fig. 7 of Martin et al. (2007) and converting IRX to 
^l(Ha), as just described, as the large, filled circles in Fig.|4] As can 
be seen, the points are in reasonable agreement with our estimates 
of the extinction as a function of stellar mass derived directly from 
Ha. If instead of using the relation we measure in Stripe 82 we take 
the relation from lCalzetti et alj feOOOl) (j3 = 0.75A(Ha) - 1.80), 
the points would be in even better agreement with t he solid line 
of Fig . [4] Using the corresponding relation fitted bv lCortese et alj 
(2006) would predict much lower Ha extinction at a given stel- 
lar mass, and be in considerably worse agreement th an the points 
shown . This may suggest that the subsample of the ICortese et alj 
(2006) survey possessing optical spectroscopy (to allow Ha, H/3 to 
be measure d) might not be representative of the larger surveys of 
iMartin et al] d2007l) (COMBO- 17) and/or the present work (Stripe 
82). However, selection bias could be present at any stage of the 
conversions we have performed (the different surveys use some 
combination of UV-selection, morphological selection, mass selec- 
tion, etc.). A detailed comparison of the IR-derived extinction is 
beyond the scope of this paper. Indeed, ideally one would like to 
have FIR data for exactly the same galaxies as all the other in- 
dicators (i.e., Stripe 82 galaxies) rather than resorting to relations 
derived from other datasets. Unfortunately, wide-field FIR imaging 
does not exist for Stripe 82 and so we cannot make a direct com- 
parison. 



6.3 Possible systematic errors 

6.3.1 Sample selection 

The main source of differences between the SFRTot-, Ha-, and 
[O II]-based SFRDs is the selection criterion. Each of the latter 
two measurements requires a S/N>2.5 in the corresponding line 
in order to be considered in the sample, whereas the SFRTot mea- 
surement is based on two different techniques (described in i]3.1| >, 
switching to the alternate technique when the emission line flux 
drops below the minimum requirement, or is otherwise corrupted. 
This is illustrated in Fig. [9] The solid lines again show the SFRDs 
based on SFRTot, Ha, and [Oil], as in previous figures (black, 
red and blue lines, respectively). The dashed black line shows the 
SFRD based on SFRTot, but now with the Ha sample's selection 
criterion (i.e., S/N(Ha)>2.5). It can be seen that this agrees very 
well with the SFRD from Ha, implying that the main difference be- 
tween the two samples is the minimum significance in the Ha line. 
The largest difference occurs around the peak of the SFRD, for 
galaxies log(M*/M0)~lO.5. This is because stellar absorption is 
strongest for the most massive (brightest) galaxies, and thus the 
lowest SFRs in Ha at these masses are hidden by underlying ab- 
sorption in the stellar continuum. Thus, it may be better to think in 
terms of an Ha equivalent width limit as well as a S/N limit. For 
these galaxies with weak SFRs and strong Balmer absorption (and 
therefore S/N<2.5), SFRTot switches to using SFRd instead of 
SFR e . From the comparison in the previous section we (and also 
ISalim et alj|2007l) have argued that SFRTot likely overestimates 
the SFR when SFRd is used instead of SFR e (i.e., for the "no- 
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Ha-SFR corrected by Balmer decrement 

SFR Tot , S/N(Hk)>2.5 

[OII]-SFR empirical correction 

- Ha-SFR, S/N([OII])>2.5 
SFRJot 




Figure 9. As Fig. [6] but imposing different S/N cuts on the sample selected. 
The solid black, red and blue lines show the nominal SFR-£ ot , Ha CO rr, 
and [O II] emp,corr SFRDs, as in Fig. [6] The dashed black line shows the 
SFRD measured from SFRTot , but now imposing the Ha selected S/N cut. 
This can be seen to agree well with the Ha-selected SFRD. The red dashed 
line shows the SFRD measured from Ha, but imposing the [O II] selection 
(S/N([0 II])>2.5). This agrees quite closely with the [0 II]-derived SFRD, 
but the latter gives a slightly higher SFRD toward higher masses. See text 
for discussion. 



Ha" or AGN component objects). The difference between the solid 
red and black dashed lines now show directly what fraction of the 
SFRD use this SFR d estimate. 

Similarly, the dashed red line shows the Ha SFRD with the 
requirement that S/N([0 II])>2.5. This brings the SFRD much 
closer to that measured by the [O II]-selected sample, albeit some- 
what lower at the high mass end. This means that the empirically- 
corrected [O II] SFRD shown in Fig. [6] will provide an excellent 
proxy for the Ha corr SFRD, provided the [O II] flux limit is suffi- 
ciently deep. 



6.3.2 Completeness, evolution and aperture correction 

We discuss in detail the effects of survey completeness and evo- 
lution in Appendix |B1 and aperture corrections in Appendix lAl For 
the present discussion, we note that the fact that our SFRD estimate 
using SFRTot in Stripe 82 is in good agreement with the SFRD 
measurement of 04 from the DR2 area means that our completeness 
and method of aperture correction give comparable results and thus 
our corrections cannot be egregious. Furthermore, our results pre- 
sented have not been corrected for evolution and our non-evolution 
corrected, 0.032 < z < 0.20 sample agrees with the evolution- 
corrected 0.005 < z < 0.22 04 sample, implying that evolution 
within our sample is a small effect. 



6.3.3 AGN contamination 

The effect of AGN contamination on the total SFRD inferred from 
emission lines such as [O II] and Ha is complicated to assess. 
04 adopted the approach of not using emission line measurements 
to measure SFR where there may be significant contribution to the 
line flux from AGN light, as determined from line ratios using the 
BPT diagram; and instead used a relation derived from star-forming 
galaxies. Furthermore, 04 use this relation for any galaxy not 
classed as star-forming (inlcuding low S/N star-forming and "no 



Ha") and apply the relation between D4000 and SFR/stellar mass 
calibrated from star-forming galaxies. Salim et al. (2007) have sug- 
gested, based on their UV-derived SFRs, that this relation used by 
04 may not be directly applicable to non-star-forming galaxies (see 
in particular their fig. 3). Since in the current work we are primarily 
interested in applying simple SFR calibrations to high redshift, we 
adopt the following. At high redshift, the common practice is to use 
whatever data is available to identify galaxies in the sample which 
may contain an AGN (e.g., MIR colours, X-ray emission; typically 
a detailed decomposition of the contribution from an AGN to the 
line flux using, for example the BPT diagram, is not possible) and 
reject them from the analysis. Thus, we use the line ratio classifi- 
cation from the current sample to reject those objects residing in 
the AGN' region of the BPT diagram and measure the SFRD in 
the same way as before. The best [O II] and Ha measurements are 
shown for this class of galaxies in the left panel of Fig.QJO] Compar- 
ing this with Fig.[6]shows that excluding AGN makes very little dif- 
ference to the total SFRD. Conversely, including AGN in the SFRD 
measurement and naively converting their line luminosities to SFR 
does not lead to a significant overestimate of the SFRD. If, however, 
only the star-forming class is considered, the total SFRD now drops 
appreciably. This is shown in the right panel of Fig. [10] It clearly 
shows that a more reasonable estimate of the total SFRD may be 
obtained by either taking all galaxies (Fig.[6j» or all but those whose 
light is dominated by AGN emission (Fig. [Toj left panel) and ap- 
plying the same SFR estimator to the resulting sample, rather than 
omitting all objects with any hint of an AGN contribution (Fig. 1101 
right panel). 

As an independent estimate of the effect of AGN, the pho- 
tometric SFR indicators (u-band and FUV) use colour selection 
to decide whether a galaxy is star-forming or not, independent 
of whether its emission lines possess sufficient S/N to assign a 
spectral class. AGN and composite (star-forming+AGN) classes 
lie predominantly in the "green valley" between the red and blue 
sequences in U V-optical CMDs o r on the red-sequence in opti- 
cal CMDs (e.g.. lSalim et"afll2007l) and thus our blue sequence se- 
lections will efficiently reject AGN. More importantly, due to the 
steep power law nature (oc v~ 15 ) of non-thermal emission from 
AGN (compared with the more Blackbody-like nature of the stel- 
lar light), the contamination to the rest-frame UV flux from AGN 
is much less problematic than the contamination to emission lines 
such as Ha and [O II], 

6.3.4 Obscured star-formation 

The main disadvantage of all the SFR indicators considered in this 
work is that a correction must be made for the effects of dust atten- 
uation and that none of the indicators is directly sensitive to dust- 
obscured star-formation. In order to directly trace dust-obscured 
star-formation, a dust-insensit ive indicator such as radio contin- 
uum emission dCondorj ll992) or Mid- or Far-IR emission (e.g., 
iKennicuttll 19981) should be employed. Unfortunately, the Stripe 82 
region does not possess sufficient overlapping radio or MIR/FIR 
data to allow a direc t test of the dust-corrections used. However, 
Hop kins et all (2003) studied a sample of 106000 galaxies in the 
SDSS Early Data Release overlapping with the FIRST 1.4GHz 
radio survey. They found that the Balmer decrement-corrected 
Ha SFR gave reasonable agreement with the radio-derived SFR, 
although the fact that only 3079 galaxies from this large sample 
had measurable radio luminosities (of which 791 were due to star- 
formation rather than AGN activity) shows the sensitivity limitation 
of using radio emission as a tracer of star-formation. 



16 D. G. Gilbank et al. 



T 




log(M./M s ) log(M,/M s ) 

Figure 10. SFRD as a function of stellar mass for different classes of galax- 
ies. Solid histogram shows SFR^i estimate, red and blue curves show 
best estimates (Balmer decrement corrected and best correction) for Ha and 
[O II] SFRs respectively. Left panel shows these measures with the AGN 
class (based on line ratios) removed, right panel shows the same for only 
those galaxies in the 'star-forming' class (after removing AGN, compos- 
ite and low S/N classes). For both panels, the same photometric measures 
of SFRD (orange dashed line shows slope-corrected FUV; magenta dot- 
ted lines shows B almer decrement-c orrected u-band; and purple histogram 
shows result from lSalim et all20 07) are shown as in Fig.[f] The photomet- 
ric samples are based only on galaxy colour to determine which galaxies are 
star-forming, and do not use the line ratio classification. See text for details. 



In i]6.2l we have shown that estimates of extinction using FIR 
results scaled from other surveys give comparable results to the op- 
tical/UV indicators we employ in this work. Obviously overlapping 
extinction-insensitive indicators such as FIR would be a desirable 
addition to the Stripe 82 data to allow a direct comparison of these 
measures. 

6.4 Application to higher redshift samples 

One of the main motivations for this work is to be able to 
use the corrections derived in the local Universe to correct the 
mass-dependence of the SFRD at higher redshift. At z~ 1, sev- 
eral surveys are now starting to produce statistical samples of 
galaxies with spectroscopic redshifts, stellar masse s and [OH]- 
derive d SFRs. The compilation presented recently in lDavies et al.l 
(2009) used three surveys to construct a sample spanning 8<J 
1or(Mi./Mq)<;12 at z~l (ROLES, GDDS, DEEP2, respec tively 
IDavies et alll2009l: I Juneau et alj|2005l : IConselice et~ai]|2007t) . The 
first two surveys, ROLES and GDDS, use only [O II] to measure 
SFR, whereas DEEP2 uses a combination of [O II] and 24/xm. For 
the moment, we just consider ROLES and GDDS' [O II]-derived 
SFR which use directly Eqn. [7j (i.e., constant luminosity conver- 
sion), but with a small correction f or IMF. For this section only 
we switch to our preferred IMF of iBaldrv & Glazebro ok (2003) 
(BG03). [O II] SFRs can be transformed from the Kroupa to BG03 
IMF using log SFRbgos = log SFRKroupa - 0.187 and stel- 
lar masses using logM t ,BG03 = log M„ y„ m „ - 0.08. This 
choice of IMF also means the IDavies et alj d2009l) numbers may 
be used directly, and we reproduce the results from their fig. 2 in 
Fig- EH The raw measurements of [O II] SFRD using the constant 
extinction prescription are shown as open circles with error bars. 



As shown in S|4] the empirical mass-dependent correction Eqn. [8] 
provides an improved estimate of the [O II]-inferred SFRD over 
the simplistic assumption of constant dust extinction. This cor- 
rection (solid symbols) lowers the SFRD in low mass galaxies 
(log(M*/MQ)<jl0) and ra ises the value in hig her mass galaxies. 
As such, the conclusion of IDavies et alj J2009) that a turnover of 
the z~l SFRD occurs below log(M*/M0)~lO is strengthened. 

Under the assumption that the SFR function has the same 
shape at z~ 1 as z~0, we can use Fig. [7j to ask if the SFRD has 
converged in the high redshift sample. In SDSS, a limiting SFR of 
~1 Mq yr~ x is sufficient at log(Af*/M0)~lO. Using our em- 
pirical correction (Eqns. [8] and |7j applied to GDDS' conserva- 
tive [O II] flux limit ( June au et alj]2005l) yields a limiting SFR of 
around 2 Mq yr~ . Similarly, ROLES would have a corrected limit 
of ~0.2 M Q jyr" 1 at log(M»/M©)~9. Fig.|7jshows that a limiting 
SFR of ~0. 1 is sufficient for the SFRD to have converged locally. 
So, both z~l surveys are quite close to this local limit. However, 
the SFRD is on average Jj3 times higher at z~l than z~0. Thus, 
if the shape of the SFR function is the same at both epochs and the 
required SFR limits are higher by a similar amount, these higher 
redshift surveys are comfortably above the SFR limit needed for 
the SFRD to have converged. 

We repeat that although we have used stellar masses from 
SDSS derived in a very different way from that commonly done at 
high redshift, if we switch to SED-fitted masses (indeed, the same 
code as used to fit the ROLES and GDDS data), our results are 
unchanged. 

It is prescient to ask if the empirical correction (Eqn. [8} 
derived at z~0 is applicable to galaxies at z~l. The correc- 
tion is likely due to the metallicity dependence of [Oil], the 
mass-dependence of the [O II]/Lto ratio and dust extinction. Re- 
cent observations of the mass-metallicity relation at z~l (e.g., 
ISavaglio et alj|2005l rCowie & Bargerll2008h suggest that the rela- 
tion is consistent with only mild evolution from that measured at 
lower redshifts, and so, since the metallicity correction is detem- 
ined via this relation, the metallicity dependence might be expected 
to be similar. Unfortunately, data to make a stringent test of all 
these factors does not currently exist. However, we can make an 
approximate check by taking the SFRD from DEEP 2 (open dia- 
monds with error bars) using 24fim SFRs dConselice et al.ll2007t) 
and comparing this with the empirically-corrected [O II] SFRs. 
Since the 24/xm emission probes radiation from star formation re- 
process ed via dust emission, th e combination of [O II]+24/im (as 
used in IConselice etahl |2007|) should be a good proxy for total 
sfrE] These can be seen to appear consistent with the corrected 
[O II] measurements in Fig. QT] Although the mass range of the 
DEEP2 data is limited to the very highest mass galaxies, this is 
where the correction should be largest and thus the agreement, al- 
beit within the broad uncertainties, is reassuring. A more thorough 
test of our correction at higher redshift will require, for example, 
overlapping Ha and [O II]-SFRs for a sample of galaxies spanning 
a range of stellar masses and such work is currently underway. 



10 In this case. lConselice et aT, ( 2007) use [O II] uncorrected for extinction 
effects, since the obscured star-formation is directly probed by the 24/im 
emission. Since the bul k of our empiric al correction accounts for dust, it 
is appropriate to use the Conse lice et alj j2007l) SFR without our empirical 
correction to [O II]. 
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Figure 11. SFRD versus stellar mass for the z~l samples from bavies et alj i20O9l) (ROLES, log(M*/M )<lO.2) and Ijuneau etal]|2005l) (GDDS, 
log(M,/MQ)> 10.2). Open circles with error bars show the [O II]-inferred SFRD from these surveys using the constant luminosity conversion [O II] SFRD. 
Filled circles show the same me asurements corrected using our empirical correction from Eqn.[8] Open diamonds with error bars show the [O II]+24/im SFRD 
from DEEP2 ( Conselice et al. 2007). The curves show the [O II]-derived local SFRD from SDSS Stripe 82 using constant luminosity conversion (dashed line) 
and [O II]corr using our empirical correction (solid line). The dotted line shows the 04 estimate for reference. 



7 CONCLUSIONS 

We have compared SFR indicators based on Ha, [O II], u-band and 
FUV luminosity using the SFRD as a function of stellar mass for 
~50000 in the SDSS Stripe 82 region (0.032 < z < 0.20). 
Our main conclusions are: 

• Using Balmer decrement-corrected Ha luminosity and a 
simple it-band based aperture correction for all galaxies irrespec- 
tive of class yields a SFRD in good agreement with the detailed 
modelling of 04. 

• In the absence of Balmer decrement information, a simple 
mass-dependent prescription (Eqn. [9} can be used for correcting 
SFR estimators for extinction. 

• [O II] luminosity works well as a proxy for Ha lumi- 
nosity, in a statistical sense, but previous empirical corrections 
underestimate the correction required toward high stellar mass 
(log(M* / Mq)> 10). We present a new empirical correction as a 
function of stellar mass which is calculated on a galaxy-by-galaxy 
basis for spectroscopically classified star-forming galaxies, and is 
thus insensitive to uncertainties in aperture correction, AGN con- 
tamination, etc.: 



SFR emP:Corr /(MQyr x ) 



L([OII])/(3.80 x 10 40 ergs^ 1 ) 



a tanh[(x — b)/c] + d 

where a = -1.424, b = 9.827, c = 0.572, and d = 1.700, for a 
Kroupa IMF. 

• Our empirical correction to [O II] is shown to reconcile the 
Ha and [O II] SFRDs as a function of stellar mass. After applying 
this correction, the [O II] SFRD is still somewhat lower (~ 30%) 



than other methods, such as Ha, but the primary reason for this dif- 
ference is the decreased sensitivity of [O II] to a given SFR thresh- 
old (set by the flux limit of the SDSS spectroscopy) at the high 
stellar mass end. 

• Photometric (u and FUV) SFRDs (for galaxies with spec- 
troscopic redshifts) give reasonable agreement with the emission 
line based SFRDs, but the effects of extinction toward high stellar 
masses become more important, particularly for the FUV. A rea- 
sonable estimate of the extinction can be made from the UV slope. 
All the photometric indicators predict a greater SFRD in low mass 
galaxies than the emission line indicators. The u-band luminosity, 
even assuming the simple constant-extinction model, is a powerful 
SFR proxy. 

• We have estimated the limiting star-formation rate thresh- 
old required for the SFRD to converge and demonstrated how this 
may be applied to higher redshift surveys, assuming that the SFR 
function in a given mass bin has the same shape as that in the local 
Universe. 

• Applying our empirical mass-dependent correction for 
[O II] SFR to an [O II]-selected sample at z~ 1 shows that the 
turnove r in the SFRD of l ow mass (log(M t /M ( v J )< > 9) galaxies 
seen by iDavies et al. I d2009h is in fact even more significant than 
reported in that paper, once the effects of metallicity and extinction 
on [O II] luminosity are taken into account. 

Mass dependent systematic effects (such as extinction) are 
of vital importance to studies of galaxy evolution. Neglecting 
such effects may lead to incorrect trends as a function of stel- 
lar mass. Using corrections derived locally, such as those pre- 
sented here, is an improvement over current methods which as- 
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sume, for example, a constant conversion between [O II] lumi- 
nosity and SFR. [O II]-selected surveys can potentially provide a 
powe rful and efficie nt probe of star-formation in the z J> 1 Universe 
(e.g.. lGilbank et alJuOlOf) . Verifying the accuracy of these empiri- 
cal calibrations will require the cross-comparison of different SFR 
indicators at higher redshifts. In particular, the inclusion of Ha, 
H/3, [OIII] and [N II] data using 8-m class, NIR, multi-object spec- 
troscopy on statistical samples will be a critical test. 
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APPENDIX A: APERTURE CORRECTIONS 

As discussed in $2.2\ it is necessary to extrapolate from flux mea- 
sured within the fibre to total flux. The total SFRs measured from 
emission lines can be dominated by the flux outside the fibre, par- 
ticularly for the most massive (and hence largest) galaxies. We have 
adopted a simple scaling based solely on the ratio of the u-band 
Petrosian flux (which should be a good approximation to the total 
flux) to the it-band flux measured within the fibre. We choose the 
it-band since this is likely to give the best proxy fo r SFR . This 
scheme is similar to that adopted by Hopkin s et al. I J2003I) . ex- 
cept that they chose to use a u-band scaling only for [O II] and 
r-band for Ha; whereas we use a it-band scaling for both mea- 
surements. Experimentation with other bands showed that using g 
or r band to aperture correct gave comparable results to using it; 
the biggest difference was that the glr aperture correction becomes 
smaller than the u-band correction by ~ 10-20% for high stellar 
mass (log(M*/M Q )^10.5) objects. 

Fig. lAll shows the aperture correction (it Petrosian flux/it fibre 
flux) as function of redshift and mass. Coloured points show galax- 
ies with stellar masses in the three mass ranges indicated. Open 
squares show the median correction for each of the three samples 
in redshift bins. It is clear that at the low redshift end of our sam- 
ple (z~0.03) the aperture corrections can be large (a factor of ~10) 
and they can be significantly larger for higher mass (larger) galax- 
ies than for the typical galaxy. By z~0.1, the median aperture cor- 
rection has converged to a value of ~3 and the median correction 
for higher mass galaxies is now not much larger than the median 
correction for galaxies of lower stellar masses, at the same redshift. 

04 derived an empirical colour-dependent aperture correction 
using [g — r) and (r — i) colours. This empirical relationship is cali- 
brated based on star-f orming galaxies, b ut applied to galaxies of all 
spectroscopic classes. [Salim et al] fc007h have cautioned that such 
a calibration does not seem to hold for non-star-forming classes, by 
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Figure Al. The aperture correction, M-band Petrosian flux/fibre flux, as a 
function of redshift. Coloured points show galaxies with stellar masses in 
the three mass ranges indicated. Squares indicate the median corrections in 
mass bins for each of the samples. These show that although the aperture 
corrections for the most massive (and hence largest) galaxies are both large 
and also even larger than the typical galaxy at lower redshifts (z<0. 1), to- 
wards the higher redshift half of our sample (0.1<z< 0.2) the corrections 
are typically a factor of three for all galaxies. See text for discussion. 



comparing 04 SFR To t with th eir GALEX-derived FUV SFRs (see 
in particular Salim et al. 2007, fig. 5). 

We wish to test whether the difference seen at the high mass 
end (Fig. [6} between our emission line measures in Stripe 82 and 
SFR.Tot from 04 (from a larger SDSS area) might be due to the 
aperture corrections used. To do this, first we directly compare 
our aperture corrections with 04's aperture corrections (by dividing 
the SFRTot value in the MPA catalogue by the SFRfn, re mea- 
surement to recover the correction used) on a galaxy-by-galaxy 
basis. We plot the log of this difference as a function of stellar 
mass in Fig. IA2I Larger filled points with error bars indicate the 
median and standard deviation of this measurement. We see that, 
on average, our aperture corrections are systematically larger by 
^0.2dex than 04, with significant scatter. There is little trend with 
stellar mass, although the scatter becomes significantly larger at 
log(M*/MQ)>10.0, with a large tail of objects towards lower dif- 
ferences (this work - 04) and at log(A:/*/Mo);>l 1.5 the median 
offset begins to increase. This trend goes in the opposite direction 
to that required to explain the higher 04 SFRD at high stellar masses 
versus ours, and so aperture corrections are unlikely to be a major 
effect. 

The fact that 04 use a colour-dependent aperture correction 
and different measures of SFR (their SFRd and SFR e ) for galax- 
ies of different spectral classes (but derived from star-forming 
galaxies) may lead to partially cancelling effects which result in us 
deriving comparable total SFRDs. We reiterate that we use a single 
aperture correction for galaxies of all spectral classes and naively 
assume that all measured line luminosity is due to star-formation. 
To check if this might be the case, we next compare SFRD esti- 
mates made entirely within the fibre (to remove the dependence on 
aperture corrections) and also examine these split by spectral type. 

Fig. IA3I shows our best estimates of Ha and [O II] SFRs 
in Stripe 82 with the aperture corrections set to unity (solid 
red and blue curves respectively). The black histogram shows 
SFRtoi measured only within the fibre (which we refer to as 
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Figure A2. Comparison of the w-band aperture correction with the colour 
dependent correction of 04. See text for discussion. 



SFRfihrc- The solid lines show these measurements just for galax- 
ies classified as star-forming. It is clear that SFR^ rc overes- 
timates the SFRD with respect to our Ha and [O II] estimates, 
even within the fibre. For the star-forming class, SFRm,i C is 
most closely based on emission line luminosity (like the Ha and 
[O II] methods) and so these should be expected to give the best 
agreement. It can be seen that even for this class, the former in- 
dicator overestimates the latter by an amount comparable to the 
difference in aperture corrections (hence the total SFRDs give 
better agreement). When all galaxy classes (rather than just star- 
forming) are considered, the agreement between the fibre-based 
measures is now somewhat better, at least for Ha compared with 
SFiiflbrc. It should be remembered that for non-star forming 
galaxies SFRabic does not directly use emission line luminos- 
ity, but rather a scaling based on SFR/stellar mass from the star- 
forming class. Out best estimate of [O II] (using the empirical cor- 
rection derived for star-forming galaxies) now underestimates the 
other two methods (as it does for total SFRDs in Fig.[6]l. Again, this 
illustrates that contamination by AGN is unimportant for [O II]- 
based measurements, and that the primary cause of the difference 
between [O II] and Ha-estimated SFRs is the sample selection de- 
scribed in £|6.3.1| 

Finally, an independent check of the validity of our aper- 
ture corrections is provided by the photometric SFR estimators. 
The agreement between our photometric SFRD estimates and our 
Ha SFRD at the high mass end reassures us that our aperture cor- 
rections are reasonable, but this is the regime where the effects of 
dust are most important and so a conspiracy between an underes- 
timate of the extinction correction in the UV coupled with an un- 
derestimate of the aperture correction in the emission line SFRDs 
might artificially cause good agreement. In higher redshift surveys 
(z~l, say), aperture corrections are relatively mild and it is only 
for this low redshift (fibre-based) comparison where the aperture 
corrections can reach factors of several for the largest galaxies that 
this is of concern. 



Figure A3. SFRD as a function of stellar mass measured within the fibre 
for star-forming galaxies. Red and blue curves show best estimates (Balmer 
decrement corrected and empirical correction) for Ha and [O II] SFRs 
and black histogram shows that using SFR^^. Only star-forming (04's 
class=l) galaxies are shown, to simplify the comparison between the in- 
dicators. The SFR measurements within the fibre are in good agreement. 
Dashed lines show the same indicators, but this time for galaxies of all 
spectral classes. The Ha and [OII]-based indicators still predict slightly 
lower SFRs than SFR^ot (which uses an alternate estimator for non-star 
forming galaxies). See text for discussion. 



APPENDIX B: COMPLETENESS IN STELLAR MASS 
AND EVOLUTION 

At a low redshift limit of z=0.005 and a magnitude limit 
of r=17.77, the SDSS is complete in stellar mass to below 
log(M*/iW0)~8 (04). For the current work, we impose a low red- 
shift limit of z=0.032 in order to include [O II] in the observed 
spectrum. The deeper Stripe 82 data impose a fainter magnitude 
limit of r=19.5. To estimate the completeness limit in stellar mass, 
we examine the stellar mass function in redshift slices. The left 
panel of Fig. IB II shows the distribution of galaxies in Stripe 82 
in stellar mass and redshift. The dashed lines indicate the redshift 
limits we impose. It is clear from the distribution of points that 
the lowest mass galaxies are only sampled at the lowest redshift 
range of the sample, due to the magnitude-limited nature of the 
sample. In the right panel, the stellar mass function is constructed. 
This follows the same Vmax procedure as the SFRD (Eqn. | lOi , 
but replacing SFR with unity to simply count objects. The log of 
the number density (log (f>) is shown as a function of stellar mass 
(\og(M* / Mq)) in Az=0.01 redshift bins. Error bars have been 
omitted for clarity. It can be seen that the overall shape of the stel- 
lar mass function over the bulk of the mass range is very similar 
for most of the redshift range. A notable exception is the very low- 
est redshift bin plotted, which shows a much lower log <j> at the 
high mass end. This can also be seen as the paucity of high mass 
objects in the left panel, and may be due to cosmic variance, due 
to the relatively low volume sampled by the lowest redshift slice. 
Nevertheless, the sh ape of the stella r mass functions are largely 
consistent (c.f., also iBell et alj 120031) and so the position of the 
turnover at low masses may be used to estimate the completeness 
limit. The turnover in the lowest redshift bin seems to occur just 
above log(M»/M0)=8.O and so the survey should be complete 
to around this limit at this redshift. This estimate may be some- 
what pessimistic, as the stellar mass function considers galaxies 
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of all spectral types/colours, whereas we primarily care about star- 
forming galaxies when estimating the SFRD. Stellar mass incom- 
pleteness will occur first in a magnitude limited survey for non-star 
forming galaxies (which have redder colours and hence higher MIL 
ratios), so this should be considered as a lower limit to the com- 
pleteness. 

The effect of evolution over the redshift range of the Stripe 82 
data is of potential concern. No correction for evolution has been 
made to our SFRDs. In order to assess this effect, Fig. IB2I shows 
the SFRD for SFRTot, and the best corrected Ha and [O II] mea- 
surements in two different redshift bins: 0.032 < z < 0.1 (left 
panel) and 0.1 < z < 0.2 (centre panel). Even if the SFRs of 
individual galaxies are not evolving, the different redshift cuts im- 
posed will result in different SFRDs, due to sampling effects, in 
the sense that the lowest mass galaxies must be drawn from the 
lowest redshift end of the survey (as they become too faint to be 
included in the magnitude selection of the survey at higher red- 
shifts). Also, since the volume probed over a fixed angular area 
survey will be larger at higher redshift, a larger fraction of the in- 
trinsically rarer higher-mass galaxies might be expected to come 
from towards the higher redshift end of the sample. The fact that 
our overall 0.032 < z < 0.20 SFRD using SFR To tis in good 
agreement with the evolution-corrected, 0.005 < z < 0.22 SFRD 
of 04 suggests that our whole sample comprises a fair representa- 
tion of the local Universe 

In the higher redshift bin, the most obvious effect is the rapid 
decline of the SFRD at log(M*/Mo)<9.5. This is obviously due 
to incompleteness, since these lowest mass galaxies are only visi- 
ble in the Stripe 82 data at the lowest redshifts in this survey. At 
higher masses, the SFRD is comparable to that seen in Fig. [6] and 
the more rapid drop off than in the 04 result is still present. In the 
lower redshift bin, the low mass tail now agrees well with that of 
Fig.|6] but the high mass SFRD is suppressed. This is likely due to 
evolution, in the sense that the universal SFR is declining towards 
the present day, and since the bulk of the SFRD occurs at higher 
stellar masses, the effect is most obvious there. 04 attempted to 
correct for evolution by assuming a universal power law evolution 
for the SFR of the form (1 + z) 3 (see 04 for the justification for 
this correction). Thus, the SFRD normalised at z—0.1, SFRD0.1, 
can be calculated as SFRD0.1 = SFRD z (l + 0.1) 3 /(1 + zf, 
for a galaxy at redshift z with an SFRD, SFRD Z . The result of 
applying this correction is shown in the right panel of Fig. IB2I It 
can be seen that this is in good agreement with the results in Fig. [6] 
and so we conclude that evolution over the redshift range discussed 
here has little impact on our conclusions. From the above discus- 
sion, assuming this form of a correction is approximately correct, 
the biggest systematic effect not applying this can have on our re- 
sults is pa 10% at the higher redshift end (and thus higher stellar 
masses), and «6% at the lower redshift end (hence lower stellar 
masses). This is comparable with the 6% effect reported by 04 (for 
a slightly wider redshift range). 
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Figure Bl. Left panel: distribution of Stripe 82 galaxies in mass and redshift. Dashed lines indicate redshift limits adopted for this survey. Right panel: stellar 
mass functions in narrow redshift slices, in order to verify the limiting stellar mass probed by our 0.032<z<0.20 sample. The turnover in the stellar mass 
function and the envelope of galaxies in mass vs redshift give an indication of the stellar mass completeness limit of the survey. Only the redshift range z<0. 1 
is plotted, since the survey clearly becomes incomplete for the lowest mass galaxies by this point. See text for discussion. 
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Figure B2. The SFRD as a function of stellar mass, as for Fig. [6] but this time split into redshift bins to study the effect of evolution. Red and blue curves 
show best estimates (Balmer decrement-corrected and best correction) for Ha and [O II] SFRs respectively. Solid histogram shows SFR^iani green dashed 
surve shows 04 SFRD determination (at z=0.1 in every case). Left panel shows 0.032 < z < 0.1, centre panel shows 0.1 < z < 0.2 and right panel shows 
the full redshift range, but SFR of each galaxy has been corrected to z=0.1, assuming SFR scales as (1 + z) 3 . See text for details. 



